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SYMBOLS LIST 


A,B,C,A^,a^,b^,c^ Constants as defined in text. 

a,b,c, RSI cell edge dimensions. 
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Eccentricity of fiber distribution 

Young's Modulus in i direction. 

Normalized number of fiber intersections per 
scanning line length per orientation of scan- 
ning line on photomicrograph . 

Continuous distribution function for fiber 
distribution in RSI. 

Discrete fiber distribution function for 
representing f(0). 

Shear modulus . 

Total length of each fiber in cell. 

Coordinate System for RSI tile 
N = through-the-thickness 
S , T = in the plane of tile. 

Least squares - sum of squared residuals. 

Variance of fiber distribution. 

Orientation of vertical planes of trusses. 

Locator of discrete angle within constant 
area interval for f(0). 

Orientation of fiber pairs in each plane. 

Thermal expansion per unit length per 
degree Kelvin. 

Fiber orientation within continuous fiber 
distribution . 

Mean orientation of fiber distribution. 

Stress component. 

Fiber orientation within discrete fiber 
distribution. 

Poisson's ratio. 

SPACE cell axes. 
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INTRODUCTION 


Analytical techniques for predicting the mechanical behavior 
of rigidized surface insulation materials (RSI) under generalized 
states of stress were developed by Materials Sciences Corporation 
(MSC) under contract NASl-11819 (the results are in Ref. 1). 
Initial computations indicated that the model was sufficiently 
realistic to be able to predict the properties of RSI material 
with reasonable accuracy, given only photomicrographs of the 
material and some basic fiber material properties. The previous 
work also demonstrated the ability of the RSI model to predict the 
sensitivity of various mechanical properties of RSI to changes in 
the numerous geometric parameters contributing to the final RSI 
material . 

The current effort has been applied to determining how a 
typical RSI material would behave under combined loading states. 

In particular, the thermal stress states induced during re-entry 
of the Space Shuttle were of prime concern. A typical RSI tile 
was analyzed for re-entry thermal stresses under computed thermal 
gradients for a model of the RSI material. The results of the 
thermal stress analyses were then used to aid in defining typical 
combined stress states for the failure analysis of RSI using the 
MSC-SPACE Program. 

An additional objective of this study was to modify the pro- 
gram described in Ref. 1 to make it more user-oriented. A com- 
plete listing of the program and supporting documentation are pro- 
vided in the Appendices to enable others to utilize the program. 


METHOD OF ANALYSIS 

The Reusable Surface Insulation (RSI) Materials studied are 
rigidized fiber networks. The analysis of RSI materials as devel- 
oped in Refs. 1 and 2 is based upon the modeling of the actual 
material by a space frame structure in which each section of fiber 
between nodes is treated as a beam having some initial curvature. 
The model utilizes a series of different, regular arrays of 
trusses, having various orientations, to represent the actual three 
dimensional fiber array, see Fig. 1. This may be regarded as re- 
presenting a continuous statistical distribution function by a 
series of discrete delta functions of differing magnitude (i.e., 
approximating a curve with a bar graph) . The analysis of RSI 
materials requires treatment of nonhomogeneous materials . When 
this material is subjected to any loading, the average deformation 
of the material depends both on the physical and geometrical vari- 
abilities of its constituents. The theory of deformation of heter- 
ogeneous media, considering physical variability from phase to 
phase has been reviewed in Ref. 2 and will not be repeated here. 
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The micromechanical RSI model follows the concept that the 
three-dimensional random network of fibers can be represented by 
a basic set of truss elements oriented in different planes, as 
introduced by Rosen and Bagchi [ 2 ] . Results obtained from the 
analysis of that model have shown encouraging agreement when com- 
pared with experimental values. This agreement has motivated the 
present development of a general model to evaluate the characteris- 
tics of RSI material exhibiting both geometrical and material vari- 
ability. The geometrical and material variabilities of the actual 
material are represented by a number of trusses oriented in differ- 
ent planes with variable material and geometrical properties. The 
concept of parallel trusses is introduced to represent such geo- 
metric factors as orientation, fiber aspect ratio, eccentricity, 
volume fraction, etc. 

Basically, the approach is to represent the continuous statis- 
tical distribution function defining fiber geometry by a series of 
elements having discrete values of each of the variables. This 
may be viewed physically as an approximate model, or it may be 
viewed mathematically as a discrete representation of a continuous 
function . 

Details of the analytical model are available in Refs. 1 and 
2. A description of how photomicrograph data of RSI material may 
be reduced to a useable form, and a description of the input vari- 
ables required for utilization of the RSI computer program, "SPACE", 
are provided in the Appendices of this report. 

The design and analysis of RSI tiles consist of several major 
aspects. Design studies require average RSI mechanical properties 
and a suitable understanding of the failure mechanisms under vari- 
ous loading conditions. Micromechanical analyses can provide aver- 
age properties based on an understanding of the microstructure of 
the material and fiber properties. In addition, micromechanical 
analyses can be utilized to develop a suitable failure criteria to 
be applied during design studies utilizing average properties and 
finite element analyses. 

Previous studies concentrated upon the mathematical develop- 
ment of the RSI model and subsequent justification of the validity 
of the model and computer program. During the current study 
several additional steps have been taken towards utilizing the 
SPACE program to develop design criteria for RSI materials. A 
finite element model of an RSI tile was utilized with average RSI 
properties to determine critical stress states existing in the RSI 
under a typical set of re-entry heating conditions . Photomicrographs 
of the RSI were analyzed to generate input for the SPACE Program 
to develop a micromechanical model of the RSI. The resulting pre- 
dicted properties of the RSI which were generated from the photo- 
micrograph data and fiber data from Ref. l, were then compared with 
experimentally measured properties and excellent agreement was ob- 
tained. The RSI model generated by the SPACE Program was then 
utilized along with the critical stress states as determined by 
the finite element analyses to develop a set of failure surfaces 
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for the RSI under combined stresses. The results, as described 
below, represent the first predictions of the behavior of RSI 
under combined stress states, based on the micromechanical model. 


RSI TILE ANALYSIS 


A typical RSI tile subjected to re-entry heating was analy- 
zed using finite elements and average mechanical properties. The 
results of the finite element analyses were subsequently used to 
define both critical stress levels and typical combined stress 
states existing within the RSI during re-entry heating conditions. 

The RSI tile modeled with finite elements is shown in Fig. 

2. The model consisted of thirty-one elements along the tile and 
twenty-four elements through-the-thickness , with symmetry boundary 
conditions applied along the centerline of the model. An aluminum 
substrate 3.81 mm thick was assumed, along with a 6.25 mm thick 
foam pad bonding the RSI material to the aluminum. The properties 
of the aluminum and foam did not significantly influence the 
analytical results because only heat-up of the tile was considered 
and the base of the tile never exceeded 75F, 

The average mechanical properties of the RSI and ceramic 
coating materials, obtained from Ref. 3, are listed in Table 1 
as functions of temperature. Through-the-thickness temperature 
gradients for a typical re-entry heat-up period were also obtained 
from Ref. 3 and are listed in Table 2. 

Several re-entry gradients were analyzed, and it was found 
that the 240 second gradient yielded the worst stress states in 
the tile. The highest stresses in the center of the tile were 
found to be the in-plane stresses several elements under the 
coating at the centerline, see Fig. 2. The stress state at that 
location was nearly unidirectional with a magnitude to approxi- 
2 

mately 0.1 MN/m in-plane, see Table 3. At the outer edge of 
the tile (Area 2 of Fig. 2) combined stresses existed with typi- 

2 

cal magnitudes of 0.014MN/m , in the through-the-thickness direc- 

■ 2 2 
tion, 0.07 MN/m in the in-plane direction, and 0.02 MN/m shear 

stresses, see also Table 3. 

The results of the finite element study of RSI under re-entry 
heatup conditions indicate that the most useful combined stress 
failure surface data would be for combined loading conditions of 
positive tension in both the in-plane and through-the-thickness 
directions, plus shear stresses. Consequently, the studies with 
the SPACE Program for combined stress failure were concentrated 
upon these types of loading conditions . 
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ANALYSIS OF RSI UNDER COMBINED STRESSES 


Before attempting to predict the behavior of RSI under com- 
bined stresses, it was necessary to define a suitable material 
model for the SPACE Program. An estimate of the fiber properties 
was available from Ref. 1, and photomicrographs of a typical RSI 
material were obtained from Ref. 4. In addition, measured mechan- 
ical properties of the RSI were available through Ref. 3, thus 
providing an experimental check on the predicted moduli and 
strengths in the in-plane and through-the- thickness directions. 
Assuming that reasonable agreement could be obtained between the 
predicted and measured properties, one could then place a reason- 
able degree of confidence in the predicted behavior under combined 
stress states. 

The photomicrographs of RSI were studied to determine fiber 
distributions in each of the principal planes of the material. 

The raw photomicrograph data consisted of number of fiber inter- 
sections per given line length versus line orientation. The fiber 
distribution function was assumed to be of the form: 

f (9) __ 1+e cos 2 (9 - 9) 

R 1 + — sin 29 

TT 

where (T gives the orientation of the distribution and e controls 
the eccentricity, or strength of orientation of the fiber distri- 
bution. The method of reducing the photomicrograph data to_a 
fiber distribution function given by equation 1 with e and 9 
known, is given in Appendix A. 

Very interesting results were obtained for the fiber distri- 
bution functions as obtained from the photomicrograph data. The 
fiber distributions for the two through-the-thickness sections, 
the N-S and N-T planes of Fig. 2, were surprisingly similar. The 
distribution functions for both planes showed 9-90° and dis- 
tribution eccentricities of 0.35 and 0.43. In contrast, the third 
plane photomicrograph, the S-T plane, showed a circular fiber 
distribution, indicating no preferential fiber orientation. 

Figs. 3 through 5 present the photomicrograph data and best 
fit distribution function, equation 1, to the data for the three 
planes examined, namely, the N-S, N-T and S-T planes . The data 
points plotted in Figs. 3 through 5 represent the data obtained 
from the photomicrographs. That is, each point represents the 
number of fibers which intersected a given length scanning line 
which was oriented at an angle of 9 + 90° at a given location on 
the photomicrograph. A fixed length scanning line was oriented 
at 9 degrees at twenty different locations on each photomicrograph 
and the number of fibers intersecting the line was determined. The 
ninety degree shift in the data was required since the highest 
orientation of fibers would show the largest number of intersec- 
tions when the scanning line was oriented at 90° to that direction. 


4 . 



Based on the above, an in-plane fiber distribution with an 
eccentricity of 0.40 and mean fiber angle of 90° was assumed 
representative of the material in the N-S and N-T planes with 
no preferential orientation in the S-T plane. Fiber truss angles 
could thus be defined as in Appendix A. Appendix B provides the 
data supplied to the SPACE Program to represent the RSI material, 
and an example of the program input and output are provided in 
Appendix C. 

A comparison of the measured mechanical properties of RSI 
and those predicted using the SPACE Program and above described 
model is provided in Table 4. As can be seen, a large variation 
occurs in the measured properties. Large experimental variations 
in properties preclude a critical comparison between the analy- 
tical model and experimentally determined properties. This does 
not prove to be a severe shortcoming for the current program, how- 
ever, because as seen in Table 3, the predicted properties fall 
well within the range of measured values. The shear modulus is 
the only property which is out of the range of the experimental 
data. This difference could arise for several reasons. First, 
the SPACE Program assumes that the fibers intersect at rigid 
joints, which would over predict the shear stiffness. Secondly, 
although test methods do exist, shear properties are difficult to 
obtain. In the case of RSI, it could be that edge effects may 
have significantly reduced the effective cross-sectional area of 
the test specimen, resulting in a substantial reduction in measured 
shear modulus . 

It is significant to note that the model accurately predicts 
the in-plane and through-the-thickness properties, which are very 
different in each direction. This result is significant in that 
the differences between the in-plane and through-the-thickness 
properties in the SPACE model have a direct correlation to the 
fiber distributions chosen to represent the spatial distribution 
of fibers. Consequently, the reasonable agreement between the 
model and measurement indicates that there is merit in the current 
method of obtaining fiber distribution data from photomicrographs 
and reducing the continuous distribution to a discrete distribution. 

Given a valid model for the RSI material, various sets of com- 
bined stress states were applied to the model to develop an under- 
standing of the behavior of the material under combined stresses. 
Although the SPACE Program employs incremental loading, the load- 
ing increments are applied to all stress components such that the 
initial stress ratios are maintained. That is, the initial stress 
tensor is increased in magnitude to failure of the material. 

Since the finite element analyses indicated that the most severe 
stress states were combinations of tension through-the-thickness 
and in-plane, along with shear, the SPACE analyses were concen- 
trated on these cases. Namely, the intersection of the failure 

surface with the a - o_ , a - a , and o - a planes for 

n s ii ns 3 ns 
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room temperature material were determined since these contours 
would reveal the basic characteristics of the failure surface. 
Figs. 6 through 8 present the intersections of the failure sur- 
face with the above-mentioned planes . Note that the regions of 
high stresses as defined by the finite element analyses and in- 
dicated in Figs. 6 through 8 fall well below the predicted 
failure surface for re-entry heat-up conditions. 


DISCUSSION 

The correlation between the analytical model for RSI and 
the measured elastic constants and strengths of the RSI in the 
two principal material directions indicates the validity of 
the analytical model. In particular, the agreement in elastic 
moduli in the two principal material directions lends credibil- 
ity to the method of defining the fiber distributions using 
photomicrograph data. 

The thermal stress analyses of the typical RSI tile provided 
some critical stress states to be evaluated using the SPACE Pro- 
gram. Since only one re-entry heating cycle was analyzed, one 
can only assume that the stress states obtained were typical and 
not necessarily "worst case". However, the intent of the thermal 
stress analyses was to provide some typical combined stress 
states, and not as an exercise to determine a "worst case" con- 
dition for a typical shuttle mission. The results of the analyses 
indicated that the most severe combined stresses were tension 
and shear in the tile. 

A large range of combined stresses did not result from the 
thermal stress analyses. Consequently, those critical stresses 
were evaluated (in-plane tension and combined tension and shear) . 
In addition, the intersection of the failure surface with the 
principal stress planes was determined to provide a better indi- 
cation of the capability of the RSI under combined stresses. 

These results indicate that a great deal of interaction between 
stress components occurs for combined stress states, requiring 
careful consideration of the failure criteria to be used for 
design purposes. 

The overall design of RSI tiles can now incorporate a reason- 
able failure criteria. That is, the SPACE Program can be utilized 
to develop a basic RSI model of a current material such that the 
SPACE Program agrees with experimental data in the principal 
material planes. One can then perform a series of computations 
using SPACE to develop an overall failure surface for the RSIr 
or the program could be utilized as a post analysis program for 
finite element results. The stress state at each finite element 
could be used as input to the program and a failure analysis per- 
formed for each RSI element. The end result would be a minimum 
margin of safety or contours of margins of safety. 
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Many modifications have been made to the SPACE Program, dur- 
ing the study, to make the program user oriented. The primary 
revisions to the program included changing the formatted input 
to namelist input, providing multiple case capability within 
each run, choosing suitable default values for many variables, 
to simplify choice of input for most RSI analyses , and improved 
documentation of the output data printed by the program. A 
description of the input to the program is provided in Appendix 
B, with a sample case input and output in Appendix C and a com- 
plete CDC 6600 listing of the program in Appendix D. 


CONCLUSIONS 


The analytical model for RSI material incorporated in the 
SPACE Program has been shown to be capable of predicting the mech- 
anical properties of RSI material with good accuracy. This achieve 
ment is significant because the RSI model is based on the micro- 
structural properties of the RSI as obtained from basic fiber prop- 
erty data and fiber geometry and distribution data obtainable from 
photomicrographs of the RSI. Thus, given some basic understanding 
of the material, one can predict mechanical properties, or con- 
versely one can use the SPACE Program to provide some guidelines 
or goals for the RSI material development. 

The combined stress failure study described earlier represents 
the first attempt to generate a failure criterion for RSI which is 
based upon the microstructure of the RSI and yet applicable to the 
design of RSI tiles from a macroscopic, average property, point of 
view . It is significant to note that no thermal stress cracking 
of the RSI studied had been observed during testing, and the com- 
bined stress study showed that, at least for the re-entry gradients 
examined, no failures of the RSI would be predicted. 

Due to continuing developments in the RSI material, the thermo 
structural design of RSI tiles should include a model for the 
failure analysis. The SPACE Program provides such a model which 
can be used as a post processor to finite element studies, conseq- 
uently, a complete CDC 6600 listing of the program and its sub- 
routines is included as an Appendix to this report. Given satis- 
factory agreement between a SPACE model and one dimensional test 
data on a given RSI material, one can simply attach that RSI model 
and the SPACE Program to a post processor for finite element analy- 
ses to develop a complete design tool for RSI tiles. 
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TABLE 1 

MATERIAL PROPERTIES 
A - RSI MATERIAL 


Temp 

K 

e n 

MN/m 2 

V e t 

MN/m 2 

g ns ^ 

MN/m 2 

V NS V NT V ST 

al/l n =al/l s =al/l t 

cm/cm/K 

116. 

40.0 

137.8 

24.1 

0.10 

-2.78E-5 

200. 

40.0 

144.7 

25.5 


-2.22E-5 

297. 

41.3 

165.4 

27.5 

1 

0 . 

478. 

46.8 

227.4 

27.5 

1 

0.61E-4 

700. 

53.7 

254.9 

27.5 

f 

1.56 

922. 

56.5 

268.7 

27.5 


2.39 

1144. 

55.8 

261.8 

27.5 


2.89 



B - COATING 

MATERIAL - 

ISOTROPIC 

Temp 

E N 2 

V NS 

al/l 

K 

MN/m 


cm/cm/K 

144. 

6.89E4 

0.20 

-0 .67E-4 

200. 

6.54 


-5.00E-5 

297. 

6.20 


0 . 

478. 

6.20 

1 

1 .06E-4 

700 . 

6.20 


2.28E-4 

922. 

4.82 

i 

3.39E-4 

1144. 

3.44E4 

T 

4.06E-4 

1422. 

6.89E2 


-2 . 22E-5 

1533. 

6.89E1 


-2.22E-5 
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TABLE 2 


TEMPERATURE DISTRIBUTION THROUGH RSI FOR VARIOUS 
RE-ENTRY HEATING TIMES 


Distance From 


Outer Surface 


Temperature, K, at 

Time 


mm 

120sec 

160sec 

2Q0sec 

240sec 

300sec 

0.00 

711 

810 

939 

1102 

1340 

2.50 

624 

696 

804 

941 

1182 

4.75 

521 

560 

621 

711 

901 

8.50 

416 

433 

456 

489 

569 

14.12 

346 

355 

365 

378 

404 

21.62 

315 

320 

325 

331 

342 

29.12 

300 

301 

302 

304 

306 

40.37 

297 

298 

298 

298 

298 

55.37 

l 

297 

297 

297 

298 

66.62 


1 

1 

I 

297 

74.12 

X 

1 

1 

1 

i 

I 

T 

OO 

297 

297 

297 

297 

297 



TABLE 

3 



SUMMARY OF MAXIMUM STRESSES FOR 
RE-ENTRY HEATING GRADIENTS IN TABLE 2 


Area of 

Stresses 

See 

Figure 2 

In-Plane 

Stress 

MN/m 2 

Through 

Thickness 

Stress 

MN/m 2 

Shear 

Stress 

MN/m 2 

Area 1 

Center of Tile 

0.10 

0 . 

0 . 

Area 2 

Edge of Tile 

0.07 

0.014 

0.02 
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TABLE 4 


MEASURED 

AND PREDICTED PROPERTIES OF 
AT ROOM TEMPERATURE 

RSI 

Property 

Measured Value 
(MN/m 2 ) 

Mean Range 

Predicted by 
SPACE (MN/m 2 ) 

In-Plane Modulus =Eg= 

=E t 227 124+275 

144. 

Through Thickness 

Modulus = E. t 
N 

46.8 26+60.6 

48.2 

Shear Modulus = GNS 

27.5 6.9+34.4 

82.7 

In-Plane Tensile 
Strength 

.200 * 

.205 

Through Thickness 
Tensile Strength 

.055 * 

.069 

Shear Strength 

.103 * 

.117 


*Design Allowable Value 
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LAYOUT OF VERTICAL TRUSSES 






TYPICAL VERTICAL TRUSS 
n = NUMBER OF "BAYS" 

= NUMBER OF SEGMENTS 
PER FIBER 
L= FIBER LENGTH 
■ LAYOUT OF HORIZONTAL TRUSS 
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Figure 1. - A typical space model. 




s 

Figure 2 - Finite element model gf RSI tile. 
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Figure 3 . 


RSI Fiber Distribution 







Figure 4. - RSI Fiber Distribution in N-S plane of the tile. 
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APPENDIX A - REDUCTION OF PHOTOMICROGRAPH DATA 


A. 1 Discrete Representation of a Continuous Distribution 

Function : 

The space frame analysis program requires representing . con- 
tinuous distribution functions with a discrete number of points. 
Spatial orientation of fibers, fiber diameters, length between 
fiber intersections , fiber eccentricities and fiber strengths 
all possess continuous distribution functions. Analytical 
representation of RSI material requires identifying a deter- 
ministic set of fiber properties which can be modeled with a 
finite number of trusses with fixed properties. The following 
section describes how a continuous distribution can be repre- 
sented by a discrete distribution of finite number of samples. 

In particular, the spatial distribution of fibers in any given 
plane is considered. However, once the continuous distribution 
function for any of the variables is defined, the discrete repre 
sentation of the variable can be determined in analogous fashion 

It is first assumed that the continuous distribution func- 
tion of a variable is known, in this case the angular orienta- 
tion of RSI fiber in a given plane. It has been assumed that 
the distribution can be represented as a first approximation 
by an elliptical distribution f(0) such that 

f{0) = - (1+e cos 20) (A-l) 

T 

where e is a measure of the eccentricity of the ellipse and 
0 is understood to be the angle of a fiber as measured from 
the mean angle or preferred fiber orientation. In practice the 
mean angle, 0, and eccentricity can be determined by using 
scanning lines, on photomicrographs of RSI material. Namely, 
the number of fiber intersections with the scanning line, 
versus orientation of the scanning line will provide an experi- 
mental distribution of f(6). A large number of scanning lines 
located at various areas of RSI material are employed to give 
a good representation of f(0), the mean value of 0 can be found 
and the eccentricity can be determined in a least squares sense, 
the results of which are described in the following section. 

The discrete representation of f(0) is given by 
N/2 , 

g (♦) = E ± <f>, (A-2) 

-N/2 W 

where N is the total number of points and = <f>^ implying 
that <}> . is normalized with respect to the mean such that the 
mean values of f(0) and g(<j>) are zero. Specific values of 
are determined as follows : 

20 . 



(A- 3) 


f (0) is divided into N regions such that 
0 , 

/ 0 f<0) de = i ; n = 1, . . . , N 

and one value of <i>^ is located within each resulting interval. 

Each <f>. could be placed anywhere within the interval 0 , to 

® n / an< ^ the general solution requires satisfying the first N+l 

moments of f(0) in order to determine each <j>^ location. However, 

a reasonable approximation of f (0) can be made by placing all 

values of i>. in the same relative position within 0 . and 0 

i c n-1 n 

and minimizing the difference between the second moments of 

f(0) and g(4>). One thus obtains the following expression to 

be solved for a, the position of the discrete variable 

between 0 , and 0 

n-1 n 


, l I { Vl + “ <«n - Vl> > 2 -s = 0 

n 2 

where S = ^ (g- ~ e) is the variance of f (0) 


(A- 4) 


9 n-l * 9 n defines the equal area intervals of f (0 ) 

N is the number of fibers = j the number of discrete 
trusses (each truss has +0 fibers) 


An example is tabulated in Table A-l giving the desired 
positions of ^ for various values of N. A graphical repre- 
sentation of two continuous distributions for e=0.3 and 0.95 
along with the discrete representations of the distributions 
are given in Figures A-l and A-2. 

The above approach is valid for a standard elliptical 
distribution. However, when one consideres how a preferred 
mean value of fiber orientation may occur during manufacture, 
a slight modification to the elliptical distribution appears 
appropriate. Assuming that the manufacturing process starts 
with a slurry of completely randomly oriented fibers, and 
then pressure is applied to compress the slurry, it is reason- 
able to assume that the fiber distribution will be symmetric 
about both the in-plane and through-the-thickness directions. 
This can be accomplished by defining f(0) such that 


f(0) = const. (1 + e cos 2(0 - je) 


(A-5 ) 
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where 


j = 1 if 0<_ 0 tt/2 

tt< 0 3ir/2 

j = -1 if ir/2 < 0 <_ it 

3tt/2 < 0 < 2 7T 

and results in the distribution shown in Figure A-3. The 
function, f(0), is normalized such that the mean value of_ 
f(0) between 0 and tt/ 2 is equal to the mean of the data, R, 

i .e. 


/ y 2 (R - f (9) ) d0 = 0 (A- 6) 

which results in 

f(6) = 1 + e cos 2 (0-j 0 ) 

1 + ~ sin 2 0 (A— V ) 

Note that, when normalized, the constant for f (0) becomes 
a function of the orientation and eccentricity of the dis- 
tribution. Although the process is slightly more involved 
algebraically, equal area intervals and discrete locations 
for placement of discrete points can be determined in the 
same manner as previously outlined. Substituting equation 
A- 7 into equation A-3 one obtains 

0. + f- sin 2 (0. - 0) = 
i 2 i ’ 

^ (1 + 2 e sin 20)-— sin 2 0 (A-8) 

TT 

where i = 1, 2, . . . N/2 

from which one can obtain 0, and define the equal area intervals. 
For most applications one ciin locate <j> . midway between 0 . and 
e i + 1' . For a 1001:6 accurate representation one can solvl for d 
in equation A-4 , substituting the variance of equation A-7 for 
S in equation A-4. In addition, data points taken from RSI 
photomicrographs_can be fit to equation A-7 in the least squares 
sense to define 0 and e; the solution involving solution of a 
pair of non-linear equations to determine 0 and e, as described 
below. 


A. 2 Least Squares Curve Fit to Distribution Functions 

Given data from the photomicrographs of RSI material which 
can be tabulated as number of fiber intersections, f., with 
scanning line versus angular orientation, 0. , of scanning line, 

one can determine the best fit values of e x and 0 in the following 
manner. 
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Define the sum of the squared residuals for equation 
A- 7 as 


N 1 + e cos 2 {0. - 0) 2 

V = Z (f. ^ } ( A "9) 

1 1 1 + 2e_ sin 2 0 

TT 


Note that here f. represents the data points divided by the 
mean value of th^ data measurements. 

Since it is desired to minimize V by proper selection of e 
and 0, we require 


3V 

3e 


- 0 


{A- 10 ) 


™ « o 
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Taking the partial derivatives, and after some manipulation, 
one can solve each equation explicitly for e in terms of 0 
as shown below: 


|V . o one obtains 
3e 


where 


A, = 


A« = 


A. « 


A„ = 


A_ ~ 


e 

= A 2 

COS 

20 + 

A^ sin 2 0 



A 3 

. 2 
sin 

20 + 

2 - 

A. cos 26 4- A- 
4 5 

sin 20 cos 

-§* 

< f i 

- 1) 

sin 

29, 

l 


-f £ 

(f i 

- 1) 

cos 

26. 

l 


-f ^ 
-i* 

sin 2 

cos 2 

20. 

i 

20. 

l 

+ 2 

TT 

"i- 1 < f i + 

1) sin 20. 

i 


(A-ll) 


tt £ sin 20. cos 20. - Z (f. +1) cos 20. 

Ill 1 


and where £ represents summation over the number of experi- 
mental data pairs, N. Note also that the f . values are divided 
by the mean value of f . as required to maintain a normalized 
distribution function. 
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From 3V 

3© 


0 one obtain 


e = ( - B + / B 2 - 4 AC) / (2A) (A-12) 

where 

A = a-^ sin 2§ + a 2 cos 29 

- 2 - 2 - 
B = b.^ sin 29 cos 29 + b 2 sin 29 + b 3 cos 29 

+ b^ cos 49 + b^ sin 49 + bg 
C = c 2 cos 29 + c 2 sin 29 


and where 

2 

a. = — Z f . cos 29. - y sin 29. cos 29. 

J. tt r i ij i i 

2 

a„ - - Z cos 29 . 

2 l 

b l * f £ “ E (f i +l) sin 29 i 

b 0 = Z f. cos 29. 

2 ii 

b, = -Z cos 29. 

^ 1 

b 4 = J ^ sin 49 ± 

b c = - n Z cos 49 . 

5 4 1 

b g = Z (f i -l) cos 29 x 

= £ E (f ^-1) sin 29^ 

c 2 = ” Z (f i -l) cos 29 i 

Solution for e and 0 can be obtained f rom_equations A-ll and A-12 
by simply substituting various values of 9 and plotting the 
resultant values of e until an intersection of the two curves 
is obtained, noting the restrictions that the derivation requires 
e 1.0, and 0 < 0 < 7 r/ 2 . 
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TABLE A-l 



5 

6 
7 
S 
9 

10 
a - 


Values of <P ^ for Representation of f (9) 
Given by Equation A-l at e = 0.95 


2 

4 

6 

8 

10 

12 

16 

20 

33.96 

8.52 

5.54 

4.19 

3.40 

2.88 

2.22 

1.93 


47.26 

22.44 

16. 32 

12.94 

10.75 

8.08 

6.60 



54.08 

30.37 

23 . 30 

19.07 

14.10 

11.36 




58.36 

35.72 

28.34 

20.47 

16.28 





61.38 

39.70 

27.41 

21.45 






63.65 

35.35 

27.01 







45.34 

33.18 







66.93 

40.35 








49.56 








70.11 

.3773 

. 3467 

.3502 

.3569 

.3639 

.3703 

.3819 

.4156 
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f (0) 


1 + 2e SIN 2 9 

ir 


j * + I if 0<9<ir/Z 
v<6< Zrr/Z 
j = “ I if ir/Z< 9<ir 
3 ^ <9<Zir 


Figure A-3 - Fiber Distribution for RSI Material 


SPACE FRAME ANALYSIS PROGRAM 
APPENDIX B 


B-I Program Description 

The SPACE Program is a space frame analysis program developed 
by Materials Sciences Corporation (MSC) for the analysis of three 
dimensional fibrous materials such as RSI. The various materials 
being developed for reusable space shuttle insulation have several 
common factors. The materials are all very low density fibrous 
systems, which, on a microscale, appear as three-dimensional arrays 
of trusses with small amounts of matrix material located at the 
fiber intersections. This observation lead to the development of 
a micromechanical model of RSI, Ref. 1,2, which can be used to 
assess the effects of various processing and geometric parameters 
on the structural properties of RSI. 

The SPACE Program treats planar arrays of fiber trusses, where 
in the planar arrays are oriented at user specified angles in space 
Fiber mechanical and strength properties are supplied by the user 
for each truss in the model, thereby allowing for inclusion of the 
effects of variable material properties in the fibers. Additional 
input data controls the geometric parameters, namely fiber length 
to diameter ratio between fiber intersections, initial fiber curva- 
tures, fiber orientations within each plane, and orientations of 
the planes in space. Finally, the applied loading, in three-dimen- 
sional stress space, is input to the program. 

Given the above data, the program computes the stress/strain 
response of the defined material under the given loading. Incre- 
mental loading is applied since the initial fiber curvature intro- 
duces non-linearity in that the effective fiber stiffness is a 
function of the applied load. 

The program prints the strain state for each load increment 
such that the nonlinear stress/strain behavior can be observed. 

In addition, each fiber pair (truss) is checked at each load incre- 
ment to determine if its ultimate stress has been exceeded. If a 
truss breaks, it is eliminated from the model, and the model stiff- 
ness and redistributed fiber stresses are recomputed at the break- 
ing load. If further trusses break', they are eliminated also, 
until either the complete model collapses, or a stable configura- 
tion is obtained. If a stable configuration is obtained, the 
loading is then increased until an additional truss fails, at which 
time the same process as above is repeated until complete failure 
of the model is attained. Details of the analytical treatment in 
the program have been described in Refs. 1 and 2 and will not be 
repeated here. A flow diagram describing the major paths of in- 
formation flow and calculations are given in Figure B-l. The re- 
mainder of Appendix B describes the input requirements of the 
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program. Appendix C gives a sample case describing the program 
output , while Appendix D provides a list of the program appli- 
cable to a CDC 6600 computer, utilizing the RUN compiler. 

Effective utilization of the SPACE Program requires deter- 
mining the constituent material geometry and properties of e 
fibers forming the RSI material. That is, the fiber proper i® s 
must be known, along with the variability of the fibers in the 
material. An effective approach to determining the spatial dis- 
tribution of fibers and representing that distribution with a 
relatively small number of trusses has been described in e ai 
in Appendix A. 

All input to the program is accomplished through a namelist 
("DATA") format for ease of keeping track of the values given to 
each variable for a given deck setup. The input can be broken 
into three main areas : control parameters determine the number 

of iterations permitted per stress solution, and define the con 
vergence and divergence criteria for the iterative solution scheme. 
Geometry definition variables comprise the major input values 
and provide a complete description of the material model and fiber 
properties. Finally, the applied loading variables define the 
initial stress state to be applied to the model, the number of 
stress increments to be permitted and the size of each stress mere' 
ment. Many of the variables have default values defined in the 
program to minimize the input required. A listing of the input 
variables along with their physical description and default values 
follows . 

It should be pointed out that the program is written such that 
multiple cases can be executed during one computer run by simply 
adding additional namelist sets with the data to be changed for 
each case. That is, the program executes the current case and 
searches the input data for an additional input namelist, if none 
exists, the program terminates. The only restriction to the above 
method of running multiple cases is that the angular orientation 
of the fibers in each truss must not be changed. That is, BT (I,J) 
may only be input in the first namelist, all other variables may 
be changed in successive runs. 


B-2 Program Input Variables 

Control Parameters 

n jiji Number of iterations permitted for determining 

strains and stresses for a particular value of 
constant applied stress. (Default value is 5) 

SEPS Convergence parameter for iterations, i.e., 

if two successive stress states differ by 
less than this value, solution has converged. 
(Default value is 0.5) 
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SUPB 


Geometry 

NOR 

NTR 

AD , BD , DD 
SGM (I , J) 


RAT ( I , J ) 
EC ( I , J ) 

EY(I,J) 

L (I ,J) 

VTI (I, J) 
TEMP (J) , 

ALPHA 

NB(I, J) 
BT ( I , J) 

AL (I) 


Divergence parameter for iterations, i.e., 
if two successive stress states differ by 
this value, solution has diverged and execu- 
tion on that load case terminates. {Default 
value is l.E+5) 

Definition Variables 

Number of vertical planes of trusses. These 
are perpendicular to one plane of fibers. 

Maximum value is four. 

Number of fiber orientations per plane. Each 
plane is assumed to have the same number of 
fiber pairs. Maximum value is thirty. 

Dimensions of the unit cell of material being 
considered. (Default values AD=BD=DD=1.0) 

Allowable fiber strength for each fiber in 

the J^ 1 truss of the 1^ plane. 

J = 1 + NTR 
1 = 1 + NOR + 1 

The (NOR +l) th plane is the horizontal plane, 

See Pig. 1. Note that the 1,2 and 3 directions 
in the SPACE Program are unrelated to the RSI 
coordinates, allowing the user to let the single 
plane of fibers, i.e., the 2-3 plane, to rep- 
resent whatever RSI orientation appears suitable. 
In the sample case, Appendix C, the major 
material distribution differences are in the N-S 
and N-T planes, consequently the model was de- 
fined so that the 1-direction was parallel to 
the N direction of the RSI, with the 2- and 3- 
directions being parallel to the S- and T- 
directions respectively. 

Fiber length to diameter ratio between joints . 

Eccentricity of fibers (initial midspan deflec- 
tion in fiber diameters) . 

Youngs modulus of fibers. 

Length of each truss (cannot exceed unit cell 
dimensions) . 

Volume fraction of each truss in the model. 

—1,3 Linear temperature variation in 1,2, 3 direc- 
tions across material element. 

Fiber coefficient of expansion, constant for 
all fibers. 

Number of bays, or intersections, per fiber. 
Angular orientation of fibers in the J^ 1 truss 

.1_T_ 

of i plane, degrees. 

Angular orientation of I^ 1 plane I=1+N0R. 
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Applied Loading Variables 


SGO (1,1) 1=1,6 Applied initial stress state on model 


H 

1! 

H 

a ll 

1=4; 

a 2 3 

CN 

li 

H 

a 22 

H 

II 

cn 
"» * 

a 13 

1=3; 

a 33 

H 

H 

a l2 


NINC 

MSCAL 


SMLT 


Maximum number of loading increments . 

(Default value is 100) 

Maximum number of re-evaluations of stress 

distributions per loading. (Used to curtail 

calculations if model cannot sustain loading 
at a given load level) (Default value is 5) 

Stress increment multiple. Stresses incre- 
ment in steps of (initial stress) * SMLT 
(Default value is 1.0) Note that SMLT in- 
creases the magnitude of the stress tensor 
and not the direction. 
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SPACE FRAME ANALYSIS PROGRAM FLOW 



Figure B-l - Space Flow Diagram 
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SPACE SAMPLE CASE 
APPENDIX C 


The sample problem describing the input data required 
to model the RSI material described eariler in this report 
is provided here, along with the output from a typical case. 

The following NAMELIST data was used to generate the 
RSI model and run two load cases: 


$DATA 

NIT=20, SEPS=1.0E~03,SUPB=2.0E03,NOR=4 ,NTR=8, 
SGM (1 , 1)=40*0.4E05, 

RAT (1.1)=40*10.0E00, 

EC(1,1)=40*1.0E00 , 

EY (1,1) =40*0. 2E0 8 , 

L(1,1)=40*1.0E00, 

VTI (1 , 1) =40* . 02E-02 , 

TEMP (1) =3*0. 0, ALPHA=0 . 0 , 

NB(1,1)=40*4.0E00, 

BT (1,1)=5*85.97, 

BT(1,2)=5*77. 81, 

BT(1,3)=5*69. 35, 

BT(1,4)=5*60.35, 

BT (1,5)=5*50. 44 , 

BT (1,6)=5*39.09, 

BT (1,7)=5*25.48, 

BT (1,8) =5*8. 99 , 

AL(1)=.OEOO, 45.0E00, -45.0E00, 90.0E00, 

SGO ( 1 , 1 ) =1 . 0,0. 0,0. 0,0. 0,0. 0,0.0, 

NINC=8 . r MXCAL=2 , SMLT=10 . , 

AD=1 . 0 , BD=1 . 0 , DD=1.0$ 

$DATA 

SGO ( 1 , 1 ) =0 . 0,0. 0,1. 0,0. 0,0. 0,0.0$ 


The following pages are an output listing of the above 
model, giving the major output data. 
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SPACE FRAME MOOtL X 


MODEL GEOMETRY 


DIMENSIONS OF MODEL A * li000006*00 8 ■ 1*000006*00 C » 1.000006*00 

VOLUME ■ 1 • 000006*00 

NUMBER OF TRUSS ORIENTATIONS ■ A 

NUMBER OF FIBER ORIENTATIONS PER TRUSS ■ B 

VOLUME FRACTION OF MQOEL - 1.249036*00 


ORIENTATIONS 

number uf 

LENGTH OF 

VOLUME 


TRUSS 

FIBERS IN 

BAYS 

TRUSS 

FRACTION 



TRUSS 





0* 

B.597006+01 

4* OOOO0E*OO 

1 * 0000VE* 00 

2,0OoOOE-04 


0* 

7*761 00E*01 

4*OOOQOE*00 

1 *000006*00 

2*000006*04 


0. 

6.935006*01 

4.000006*00 

1 .000006*00 

2,000006-04 


0. 

8.035006*01 

4;oooooe*oo 

1.000006*00 

2,000006-04 


0* 

5,044006*01 

4,00000£*00 

1.000006*00 

2.000006-04 

O to 

0. 

3, 90900E+01 

4i OOOOOE* 0 0 

1 , OOOOOE* 0 0 

2.00000E-04 

to 

0. 

2,549006*01 

4,OOOOOE*CO 

1 , 000006*00 

2 . 000006-04 

to to 

0. 

8.990006*00 

4.00000E*00 

i*oooooe*oo 

2,000006-04 

£ to 

4.5oO00E*01 

8.597006*01 

4, 000006*00 

1 * 000006*00 

2.00000E-Q4 

to 

4.500006*01 

7,781006*01 

4.00000E*00 

1,OOOOOE*00 

2.000006-04 

to s 

4.500006*01 

6.935006*01 

4,oooooe*oo 

1,000006*00 

2.00000E-04 

c 2 

*■5000 OE *01 

6, 035Q06 ♦ 0 1 

4, 00 0 0 OE*00 

1.00000E*00 

2, OOOOOE-04 

to S3 

*.50000E*01 

5, 0*400t *0 1 

4.000006*00 

i,oooooe*oo 

2, OOOOOE-04 

to G 

4.50000E*01 

3*9090OE*0l 

4.000006*00 

1*000006*00 

2* OOOOOE-Q4 

£ £ 

4.50000E*01 

2.548006*01 

4, 000006*00 

1 . U0QQOE*O0 

2* OOOOOE*C4 

X to 

4.500006*01 

8.990006*00 

4,oooooe*oo 

l,0000OE*0O 

2*000006-04 

£? to 

-4.50000E*01 

8.597006*01 

4, OOOOOE*flO 

1.000006*00 

2. OOOOOE-04 

to to 

-4.500006*01 

7,761 00t *0 1 

4, OOOOOE* 00 

1.000006*00 

2. OOOOOE-04 

to k? 

•4.500006*0 1 

6.935006*01 

4.000006*00 

1,000006*00 

2,000006-04 

CO ^ 

•4,500006*01 

6.035006*01 

4, OOOOOE* 00 

1 , OOOOOE* 00 

2*000006-04 

to O 

-4,500006*01 

5, 04400E*0 1 

4, OOOOOE* 00 

l.OOOOUE+OO 

2. 000006-04 

O to 

-4,500006*01 

3,90900t*Ol 

4,000(M>E*00 

1.000006*00 

2.000006-04 

to 

-4.500006*01 

2.548006*01 

4, OOOOOE 4 00 

1,000006*00 

2*000006-04 

£ to 

-4.500006*01 

8,990006*00 

4, OOOOOE* 00 

1 ,000006*00 

2.000006-04 

rO to 

9* OOOOOE* 0 1 

8,59700t*0l 

4, OOOOOE* 00 

1.000006*00 

2. 000006-04 

to 

9 * OOOOOE* 0 1 

7*761 006*0 1 

4* OOOOOE* 00 

1 • O00Q0E*00 

2*000006-04 


9.000006*01 

6.935006*01 

4, OOOOOE* 00 

1,000006*00 

2,000006-04 


9,000006*01 

6.035006*01 

4,000O0E*00 

1.000006*00 

2.000006-04 


9» 000006*01 

5,044006*01 

4, OOOOOE* 00 

1.000006*00 

2. OOOOOE-04 


9.000006*01 

3, 909006* 0 1 

4* OOOOOE* 00 

1.000006*00 

2,000006*04 


9*000006*01 

Z, 54»<J0E* 01 

4, 000006*00 

1 * 000006*00 

2.000006*04 


9,000006*01 

8,990006*00 

4,00000£*00 

1.0000U£*00 

2* OOoOOE-04 



8*597006*01 

4*000006*00 

1.000006*00 

2.00000F-04 



7,781 006*01 

4, OOOOOE* 00 

1 ,0000^6*00 

2. OOOOOE-04 



6.935006*01 

4, OOOOOE* 00 

1,000006*00 

2 . OOOOOE-04 



6. 036006*01 

4, OOOOOE* 00 

1 * 000006*00 

2.000006-04 



5, 044006* Ol 

4,0OOOOE*OO 

1,000006*00 

2.000006-04 



3.909006*01 

4.000006*00 

1*000006*00 

2.00000F-04 



2,54*006*01 

4, OOOOOE* 00 

1.000006*00 

2.000006-04 



B, 990006*00 

4.000006*00 

1.000006*00 

2.000006-04 



FIBtR PROPERTIES 




ORIENTATIONS 

YOUNG I S 

ECCtNI . 

CROSS- 

MAXIMUM 


TRUSS 

FIRFRS in 

MODULUS 


SECTIONAL 

SIRESS 



TRUSS 



AREA 



uj 0 * 

8.597006*01 

2, OOOOOE* 07 

1 « OOOOUE+ 00 

2 » 4B463E-02 

4,000006*04 


<* 0* 

7,761 006*01 

2, OOOOOE+OT 

1 , OOOOOE'* 00 

2.75239E-03 

4,000006*04 


0. 

6.935006*01 

2.00000E+07 

1, 000006*00 

9.06737E-O4 

4,000006+04 


0. 

6.035006*01 

2 • 00000E+07 

1,000006*00 

5, 0 1 4476-04 

4,000006+04 


0* 

5.04*006+01 

2.00000E+07 

1.0000U6+00 

3.025436*04 

4, OOOOOE* 04 


0, 

3.909006*01 

2 « 000006* 07 

1.000006*00 

2.037106-04 

4,000006*04 


0. 

2 • 546 o Ot * 0 1 

2.000006*07 

1 , OOOOOE* 00 

1.505B7E-04 

4,000006*04 


0 ■ 

8,99000t*00 

2* 0 00 0 0E4 07 

1 , OOOOOE 4 00 

1.25T90E-04 

4.000006*04 


4.50000E+01 

8, 597006+01 

2.00000E+07 

1 ,000006 + 00 

2.464836*02 

4,000006+04 


a.sooooe+oi 

7.781006*01 

2.00000E+07 

l,O0000E*0O 

2.75239E-03 

4, OOOOOE* 04 


4,500006+01 

6,935006+01 

2. 000006+07 

1,000006*00 

9.8673TE-04 

4,000006+04 


4,500006*01 

6, 035006* 0 1 

2,000006*07 

1.000006*00 

5,014476-04 

4,000006+04 


4,5ooCOE+01 

5.044006*01 

2.00000E*07 

1.000006*00 

3.025436-04 

4,000006*04 


4.50000E+01 

3.90900E+01 

2.00000E+07 

1 , 00QO0E*0O 

2.037106-04 

4,000006*04 


4.50000E+01 

2.54R006+01 

2.00000E+07 

l.oooouc+oo 

1 . 5Q5B7E-04 

4, 000O0E*04 


4.50000E* 01 

B.99000E+00 

2 ^OOOOOE* 07 

1,000001*00 

1.2579GE-04 

4,000006*04 


-4.500006*01 

8 .597006+01 

2.000006*07 

l .OOOOOE* 00 

2.464636-02 

4.000006*04 


•4,500006*01 

7 ,781 006+01 

2 ,OOOOOE4 07 

1.000006*00 

2,75?39E-03 

4.000006*04 


-4,500006*01 

6.935006+01 

2, OOOOOE+07 

1.000006*00 

9,867376-04 

4,000006*04 


-4.500006*01 

6,035006+01 

2 • 00000£*07 

1,000006+00 

5. 014476-04 

4,000006+04 


•4 .50COOE+ 0 1 

5.044006+01 

2,000006*07 

1,OOOOU£*00 

3,02543E*04 

4, OOOOOE+04 


-4.5OO00E+01 

3,909006*01 

2.0000QE+07 

l.OOOOOE+OO 

2.03710E-04 

4, 000006+04 


-4.50000E+01 

?,54AoO£+OI 

2. OOOOOE* 07 

1 • OOOOOE* 00 

1 . 505676-04 

4,000006*04 


-4.5OOOOE+01 

8,990006+00 

2*000006+07 

1,000006*00 

1,257906-04 

4.000006*04 


9*000006*01 

8,597006*01 

2 • 000006+07 

l i 000006*00 

2.48463E-02 

4,000006+04 

*9 in 

9;OOOOOE*01 

7 • 781006*01 

2 ,000006+ 07 

1.000006*00 

2,752396-03 

4,000006+04 


9.000006*01 

6. 935006+ 01 

2, OOOOOE* or 

l.OOOOVE+OO 

9. 867376-04 

4. OOOOOE* 04 

C? « 

9.000006+01 

6.035006*01 

2. OOOOOE+07 

1 . 00000E*O0 

5.01447E-04 

4,000006+04 

& to 

9.00000E+01 

5, 044006*01 

2, OOOOOE* 07 

1.00000E*00 

3.02543E-04 

4*000006+04 

£ o 

9. OGOOOE + Ol 

3.90900E+01 

2 . 00000E*07 

1 .OOOOOE* 00 

2.037106-04 

4, O0000E+04 

it* 

9.000006*01 

2,548006*01 

2,00000t*07 

1 , OOOOOE* 00 

1,505876-04 

4. 000006+04 

& £ 

9 .OOOOOE* 01 

8*990006+00 

2*000006+07 

1.000006*00 

1.25790E-04 

4.000006+04 

JJ7 


8.597o06*0l 

2*000006+07 

1.000006*00 

Z. 484636-02 

4.000006*04 

Ay /Q 


7,781 O06+Q l 

2* OOOGOE + 07 

1.000006*00 

2,752396-03 

4,O00O0E*O4 

A^ t-7 


6,935006*01 

2, OOOOOE*07 

1.000006*00 

9.B6737E-04 

4, OOOOOE* 04 

,N 


6.035006*01 

2,OOOOOE*OT 

l.OOOOOE+OO 

5.01447E-04 

4,000006*04 

£*7 '■-. 


5,044006*01 

2, OOOQOE*07 

l.OOOOUE+OO 

3.02543E-Q4 

4.000006+04 

> '-) 


3.909006*01 

2, OOOOOE* 07 

l.OOOOOE+OO 

2.0371OE-O4 

4.000006*04 

£ N 


2.548006+01 

2,0O00OE*O7 

1,000006*00 

1,505876-04 

4,000006+04 

*© O 


8.990006*00 

2, OOOOOE* 07 

l.OOOOUE+OO 

I.25790E-O4 

4,000006+04 

£5 ^7 









ANALYSIS 





SOLUTION FOR STRESS CONVERGES IN 2 

ITERATIONS 







DIRECTION 

APPLIED LOAD 

STRAIN 



11 

1.000006+00 

1,41636E*04 



22 

0. 

-3. 81 2636*05 



33 

0. 

-1.36260E-05 



12 

0. 

-3,n030E-13 



23 

0, 

3.909766-19 



31 

0, 

-1.34076E-1 8 



orientations 

STRESS 

OWItNTATIONS 

STRFL' 

TRUSS 

FIREPS IN 

TRUSS 

FIBERS IN 



TRUSS TRUSS 



'LZ 



-3.69851E*02 ! 

1 .63463L+02 ! 

l. 4l896t* 03 1 

3.24466E*Q3 1 

-3.826196*02 2 

*1 *58043fc*02 2 

3.708106*02 2 

1 » 14o2{Jt*o3 2 

-3.82619t*0Z 3 

-1.5804Jfc*02 3 

3.708l0E*02 3 

1#14020E*03 3 

-3.953076*02 4 

~4 # 79S5QE«02 4 

-6.777526*02 4 

-9'66187t*0Z 4 

-1.07l07fc*03 5 

-4,536086*02 5 

1 ,000336*03 8 

3.114446*03 5 


SOLUTION FOR STRESS CONVERGES IN 2 ITERATIONS 


direction 

u 

22 

33 

12 

23 

31 

ORIENTATIONS 
TRUSS FIBERS IN 

TRUSS 


APPLIED LOAD 

1 • 100006*01 
0 - 
0 * 

0 . 

0 . 

0 . 

STRESS 

TRUSS 


1 1 

1 3 

1 5 

1 7 

2 I 

2 3 

2 5 

2 7 

3 1 

3 3 

3 5 

3 7 

4 1 

4 3 

4 5 

4 7 

5 1 

5 3 

S 5 

5 7 

MEMBER FAILURE I TRUSS ] 

MEMBER FAILURE t TRUSS 5 


•4.05825E*03 
1.779956*03 
1 #547366*04 
3.5238+6*04 
-4.198l2t*03 
-1.742076*03 
4.038956*03 
1*242386*04 
•4* 190126*03 
-1*742076*03 
4.038956*03 
1,2423BE*04 
•4*337996*03 
-5.264156*03 
-7*445046*03 
•1*061836*04 
-1.17732E*04 
•5.011106*03 
1*009276*04 
3*362696*04 

F106R ft 
FIBER A 


1 

1 

1 
1 
? 

2 
2 
2 
3 
3 
3 

3 

4 
4 
4 

4 

5 
5 
5 
5 


SOLUTION FOR STRESS CONVERGES IN 2 ITERATIONS 


2 

4 

6 

8 

2 

4 

6 

8 

2 

4 

6 

8 

2 

4 

6 

0 

2 

4 

6 

6 


-1.920056*02 
7 .000346*02 
2.296966*03 
3.960796*03 
-*3.080666*02 
6.828726*01 
7.40771E*02 
1**42116*03 
-3.080666*02 
6.828T2E*01 
7.407716*02 
1.442116*03 
-4.233276*02 
-5.643696*02 
-8* 16432E*02 
-1.079406*03 
-8.660916*02 
1 . 686066*02 
2.01707E*03 
3.94362E*03 


STRAIN 


1.55050E-03 
•4. 19227E-04 
-1.49514E-04 
-3.39379E-12 
5*739416*18 
-1*06490C*17 

ORIENTATIONS STRESS 

FltfENS IN 
THU55 


Z 

4 

6 

8 

2 

4 

6 

8 

2 

4 

6 

0 

2 

4 

6 

a 

2 

4 

6 

e 


-2.120096*03 

7.651346*03 

2.500046*04 

4.29440E«04 

-3.302836*03 

7.333976*02 

8.074736*03 

1.570546*04 

-3.30283E*O3 

7.333976*02 

8.074736*03 

1.570546*04 

-4.645466*03 

-6.197486*03 

-8.970816*03 

-1.186356*04 

-9.520746*03 

1.806046*03 

2.195176*04 

4.275866*04 




o 

r; ^ 


a* 

L*3 



zmw 


Ul 

00 


h} 

M 



DIRECTION 

APPLIED LOAD 


STRAIN 



11 

1 • Q2448E*0l 


2»38060E"03 



22 

0* 


-6.38087E-04 



33 

0* 


»2*27040€«04 



12 

0* 


-5*2D077E-12 



23 

0* 


7,3SilOE-18 



31 

0* 


-1.96920EM7 



ORIENTATIONS 

STRESS 


ORIENTATIONS 

STRESS 

TRUSS 

FIBERS IN 

TRUSS 


FIBERS IN 



TRUSS 



TRUSS 


1 

1 

«6.15B39E*03 

1 

2 

•3* 16578E 4 03 

1 

3 

2.79563£*03 

1 

4 

1 • 179Q3E* 04 

1 

5 

2.37521E*C4 

1 

6 

3* 02882E* 04 

1 

7 

5,307Q4E*O4 

1 

8 

0. 

2 

1 

-6.37Z71E*C3 

? 

2 

-5.12078£*03 

2 

3 

-2.60118E*03 

? 

4 

1 • 2OO36E+03 

2 

5 

6»27354E*03 

2 

b 

1.24615E+04 

2 

7 

1 «91 224E+04 

? 

8 

2#41435E*04 

3 

1 

-6.37Z71E*C3 

3 

2 

•5* 12078E*03 

3 

3 

-2 f 60U8E*C3 

3 

4 

1.20036E*Q3 

3 

5 

6.Z7354E*o3 

3 

6 

1*24615E*04 

3 

7 

1.91Z24t*04 

3 

8 

2.4143SE*04 

4 

1 

-6,58703E*03 

4 

2 

-7 * 05551 £*03 

4 

3 

*7.99820E*03 

4 

4 

-9.42025E*03 

4 

5 

M.132l0t*04 

4 

6 

-1#36455E*04 

4 

7 

-1.61552E*04 

4 

8 

»1 • B05ZI E*04 

5 

1 

-1.79121E*04 

■5 

2 

- 1 • 44737£*04 

5 

3 

-7.55214t*03 

5 

4 

2# B9387E*03 

5 

5 

1.67954E404 

S 

6 

3.3670SE*04 

5 

7 

5« 1 73B5E* 04 

5 

8 

0* 

MEMBER 

FAILURE! TRUSS 1 

FIBER 7 





SOLUTION for STRESS CONVERGES IN 2 ITERATIONS 


DIRECTION 


APPLIED LOAD 


STRAIN 


11 

22 

33 

12 
23 
31 

ORIENTATIONS 

TRUSS FIBERS IN 

TRUSS 


1.02448E*0l 

0* 

0» 

0* 

0* 

0. 

STRESS 


TRUSS 


3.04913E-03 

•6.27776E-04 

-2*59670E-04 

-6.61388E-12 

9,12515€-ia 

-1.96952EM7 

ORIENTATIONS STRESS 

FIBERS IN 
TRUSS 


I 

I 

1 
1 

2 
2 
2 
2 
3 
3 
3 


1 

3 

5 

7 

1 

3 

5 

7 

1 

3 

5 

7 


-6.996B9E+03 t 

4* 36246E *03 1 

3* 08969E* 04 I 

0, 1 

-7* 2721 3E* 03 ? 

-2.56656fc*03 ? 

8,5C065t*03 2 

? * 4502 It * 04 ? 

-7 * 27213E * 03 3 

-2*56656E+03 3 

R« 50065E* 03 3 

3 . A t,n? i ► * r>4 * 


2 

4 

6 

8 

2 

4 

6 

8 

2 

4 

b 

4 


-3*2251 OE*03 
1.57600E*04 
4#92602E*04 
0* 

-5*71 0 18E* 03 
2 • 1 7590£* 03 
1 • 02099E* 04 
3»07484F*Q4 
-5» 71 0 18E*03 
?• 17590E*O3 
U8?099E*04 

-a rtT <L a A C 4 n 4 l 



4 

l 

-7.54736t*l>3 

4 

3 

*9* 49767E* 03 

4 

5 

- 1 * 40896t* 04 

4 

7 

»2, 0769SE* 04 

5 

1 

•2,32309t* 04 

5 

3 

•9»9371 3E* 03 

S 

5 

2,13090E*04 

5 

7 

6.5902O£*O4 

MEMBER FAILURE* 

TRUSS 5 

FIBER 7 


SOLUTION FOR STRESS CONVERGES IN 2 ITERATIONS 


DIRECTION 


11 
22 
33 

12 
23 
31 

ORIENTATIONS 

FIBERS IN 
TRUSS 


APPLIED LOAD 


l*02448E 4 01 

0* 

0 . 

0. 

0. 

0. 

STRESS 


1 1 

1 3 

1 5 

1 7 

2 1 

2 3 

2 S 

2 T 

3 1 

3 3 

3 5 

3 7 

4 1 

4 3 

4 5 

4 7 

5 l 

5 3 

5 5 

5 7 

MEMBER FAILURE t TRUSS 1 


-9.79426E*03 
5*S8877t *03 
4,1 424 lE*04 
0* 

•1.016OUE*O4 
•3, 61865E* 03 
1,1 756BE+ 04 
3,39367E*04 
M.0160I>E>04 
-3,61B65E*03 
1 « 17568E+ 04 
3* 39367E* 04 
-l,05257t^04 
•1 ,28291E*04 
-1.82519E*04 
-2,6I399E*04 
-2.89745E*04 
•1 , 1 3004E+ 04 
3* 01517E+ 04 
0. 

FIBER A 


u> 

k£> 





• 8* 19484E* 03 
• 1 ♦ 146Z9E*04 
-1.73017E*04 
-2,33904E*04 
•l.H8245E*04 
3,477S7E*03 
4* 29090E 4 04 
0. 


STRAIN 


TRUSS 


4.UBT5E-03 
-1,033B?E-03 
-3,52777E-<}4 
-8,99540E-12 
1 » 273B5EM7 
-2.23711E-17 

OHIENTATION5 STRESS 

FIOEH5 IN 
TRUSS 



-4.68641E*03 
2*1001 3E 4 04 
6.61308E*04 
0 • 

«.7*98879E*03 
2* 97420E 4 03 
2,24521E*04 

4 »25?60E 4 04 
-T,90B79E*Q3 
2*97420E 4 03 
2. 24521 E*04 
4,25760E*04 
-1.12904E*0^ 

• 1 ,5150oE 4 04 
-2.2045lE*04 
*»2*92347E # 04 

• 2* 31 0982*04 
6* 52201E* 03 
5*8feB75E*04 
0* 



SPACE FRAME MODEL 2 


MOUEL GEOMETRY 


OIMENSIONS OF MODEL A * i • 000006+00 6 ■ 1 * 000006+00 C “ i. 000006+00 

volume * i,oooooe*oo 


NUMRER OF TRUSS ORIENTATIONS « 4 

NUMBER OF FIBER ORIENTATIONS PER TRUSS = 0 

VOLUME FRACTION OF MODEL ■ 1.24903E+00 




ORIENTATIONS 

NUMBER OF 

LENGTH of 

VOLUME 



TRUSS 

F1RER5 IN 

BAYS 

TRUSS 

FRACTION 




TRUSS 






0 * 

6,597006*01 

4.000006+00 

1 « OOOOOE + OO 

2. 000006*04 



0 * 

7.761006*01 

4. OOOOOE+OO 

l.OOOOUE+OO 

2. 000006-04 



0 * 

6.93500E*01 

4, OOOOOE+OO 

1. OOOOOE+OO 

2.000006*04 



0 . 

6. 03S00t*0 1 

4 * OOOOOE + OO 

1, OOOOOE + OO 

2.00000E-Q4 



0 . 

5. Q4400E+0 1 

4,000006+00 

l.OOOOOE+OO 

2.QQOOOE-04 



0 . 

3.90900E+01 

4,000006*00 

1 . OOOOOE+OO 

2.000006-04 



0 . 

2.54800E+0 1 

4.000006+00 

l.OOOOOE+OO 

2.000006-04 



0 . 

8*990006*00 

4,000006+00 

l.OOOOOE+OO 

2.000006-04 



4.500006+01 

8.59700E*01 

4,000006+00 

l.OOOOOE+OO 

2.00000E-04 



4 • 50000E + 0 1 

7*78] 006*01 

4,000006+00 

l.OOOOOE+OO 

2. OOOOOE-04 



♦ .50000E* 01 

6.93800E+01 

4.00000E+00 

1, OOOOOE+OO 

2.000006-04 



4.50000E+ 01 

6,O350OE*Ol 

4,000006+00 

l.OOOOOE+OO 

2. OOOOOE-04 



4* 50000E+ 01 

5.044006*01 

4. OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 



4.5oOOOE*Ol 

3.909006*01 

4,000006+00 

l, 000006+00 

2*000006-04 



4 ,500006+ 0 ) 

2.540006*01 

4, OOOOOE+OO 

l.OOOOOE+OO 

2. OOOOOE-04 



4.500006+01 

6,990006*00 

4, OOOOOE+OO 

l.OOOOOE+OO 

2, OOOOOE-04 

c 


-4 ,50000E+0 l 

8,597006*01 

4. OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 

jS 

-4,50000E*01 

7. f8l006*0l 

4, OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 

A 

-4.500006+01 

6,93500E*01 

4, OOOOOE+OO 

1 .OOOOOE+OO 

2. 000Q0E-04 


-4.50000E+0 1 

6.035006*01 

4, OOOOOE* 00 

l.OOOOOE+OO 

2. 000006-04 


•4 .50000E*0 1 

5.044006*01 

4.000006+00 

1 ,000006+00 

2. 000006-04 

A* t 


•4 .50000E+0 1 

3.90900E+01 

4. OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 

<N A 

J 

-4.500006*01 

2,54 &00E* 0 1 

4. OOOOOE+OO 

1, OOOOOE+OO 

2,000006-04 


-4.500006*01 

8*99000E*00 

4, OOOOOE+OO 

1 , 00000£*00 

2 ,000006-04 



9.00000E*01 

0.597006*01 

4, 000006+00 

l.OOOOOE+OO 

2. 000006-04 


9 • 000006*01 

7.781006*01 

4, OOOOOE+OO 

] .OOOOOE+OO 

2, OOOOOE-04 


9 .OOOOOE * 01 

G.93500E*01 

4. OOOOOE+OO 

1 ,000006+00 

2,000006-04 


9*OOOOOE *01 

6.035006*0) 

4,000006+00 

1 .OOOOOE+00 

2 .000006-04 



9* OOOOOE*0 1 

5.04*OOE*01 

4, OOOOOE+OO 

1 .000006+00 

2.000006-04 


9.000006*01 

3.909006*01 

4. OOOOOE+OO 

1. OOOOOE+OO 

2, OOOOOE-04 

<9 * 


9.000006*01 

2.54fl00E*01 

4,000006*00 

l.OOOOOE+OO 

2,000006-04 


9.00000E*01 

0,990006*00 

4, OOOOOE+OO 

l.OOOOOE+OO 

2*000006-04 



8.59700E*01 

4,000006+00 

l.OOOOUE+OO 

2.000006-04 

ft/ ^ 



7.701006*01 

4. OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 




6.93500E*01 

4,000006+00 

1. OOOOOE+00 

2, 000006-04 




6.035006*01 

4, OOOOOE+OO 

l.OOOOOE+OO 

2,000006-04 



5.044006*01 

4.000006*00 

] .000006+00 

2.000006-04 




3,909006*01 

4. OOOOOE+OO 

l.OOOOOE+OO 

2.000006-04 




2,545006*01 

4, OOOOOE+OO 

l.OOOOUE+OO 

2,000006-04 




8,990006+00 

4. OOOOOE* Of) 

1 .UOOOOE + OO 

2,000006-04 


FlfitK PROPERTIES 



ORIENTATIONS 

ydungis 

ECCtNT* 

CROSS* 

TRUSS 

FIBERS in 

MODULUS 


sectional 


TRUSS 



AREA 

0* 

B.59700&+01 

2.000006*07 

1.000006*00 

2,484636-02 

0* 

7*701006*01 

2*000006+07 

i*oooooe*oo 

2-752396-03 

0* 

6*935006+01 

2,000006+07 

1*000006*00 

9,867376-04 

0* 

6*035006+01 

2,000006+07 

l.OOOOOE+OO 

5.014476-04 

o. 

5*044006+01 

2,000006+07 

1,000006+00 

3-025436-04 

0* 

3*909006+01 

2,000006+07 

1 * OOOOOE + OO 

2*037106-04 

0* 

2*54«006+0l 

2,000006+07 

l.OOOOOE+OO 

1-505876-04 

0 , 

8*99O0Ot+00 

2,000006+07 

1 * OOOOOE + OO 

1.257906-04 

4, 500006*01 

6. 597006+01 

2,000006+07 

1.000006+00 

2.484636-02 

4*500006+01 

7.761006+01 

2,000006+07 

1*000006+00 

2.752396-03 

4*500006+01 

6,935006+01 

2,000006+07 

1*000006+00 

9-867376-04 

4-500006+01 

6*035006+01 

2,OOQOOE*07 

1,000006+00 

5-014476-04 

4-500006+01 

5*044006+01 

2,000006*07 

i. oooooe+oo 

3.0Z543E-04 

4.500006+01 

3,909006+01 

2.000006+07 

1.000006+00 

2.037106-04 

4,500006+01 

2.548006+01 

2,000006*07 

1.000006+00 

1-505876-04 

4,500006+01 

6,990006+00 

2*000006+07 

1*000006+00 

1-257906-04 

-4.50000E+01 

8*597006+01 

2.000006+07 

1 * 000006 + 00 

2*484636-02 

•4*500006+01 

7. 781 006+01 

2,000006+07 

1 . Q00006+00 

2- 752396-03 

•4.50000E+01 

6.935006+01 

2-000006+07 

1,000006+00 

9.86737E-04 

•4-500006*01 

6*035006+01 

2*000006+07 

1,000006+00 

5-014476-04 

•4.50000E+0 1 

5.044006+01 

2,000006+07 

1.000006+00 

3,025436-04 

•4*500006+01 

3.909006*01 

2*000006+07 

1.000006+00 

2.037106-04 

-4.500006+01 

2*546006+01 

2-000006+07 

1*000006+00 

1.50587E-04 

•4.50000E+01 

0. 990006+00 

2.000006+07 

l.OOOOOE+OO 

1.257906-04 

9.00000E+01 

6,597006+01 

2-00000E+07 

1,000006+00 

2.4B463E-02 

9.000006+01 

7.781006+01 

2*000006+07 

1 .OQOOOE+OO 

2.752396-03 

9-OOOOOE+Ol 

6*935006+01 

2- OOOOOE+07 

1 *000006+00 

9. 867376-04 

9.00000E+01 

6*035006+01 

2.00000E+07 

l.OOOOOE+OO 

5.01447E-Q4 

9*000006+01 

5.044006+01 

2.000006+07 

1*000006+00 

3.02543E-O4 

9*000006+01 

3.909006+01 

2 .00000E+07 

l.OOOOOE+OO 

2.03710E-04 

9 *000006+01 

2.54R006+01 

2.0O000E+07 

l.OOOOOE+OO 

1.5O507E-O4 

9*000006+01 

8.990006+00 

2,000006+07 

l.OOOOOE+OO 

1.25790E-04 


8.597006+01 

2,000006+07 

l.OOOOOE+OO 

2. 48463E-02 


7*781006+01 

2.000006*07 

l.OOOOOE+OO 

2.75239E-03 


6*930OOE+O1 

2.000006+07 

1.000006+00 

9.867376-04 


6.035006+01 

2.00000E+07 

1*000006+00 

5.014476-04 


5,044006*01 

2. OO000E + 07 

1 *000006+00 

3-02543E-04 


3.90900E+01 

2. OOOOOE + 07 

1.000006+00 

2. 037106-04 


2.54Ao0t+0l 

2;oooooe+o7 

1,000006+00 

1.505876-04 


B.99000E+00 

2,000006+07 

1. OOOOOE+OO 

1.257906-04 


ANALYSIS 



solution for stress converges IN 2 

iterations 






DIRECTION 

APPLltn LOAD 

STRAIN 


11 

0- 

-1.362946-05 


22 

0- 

-9.37055E-O6 


33 

1-000006+00 

4.78628E-05 


12 

0- 

8.52167E-14 


?3 

0. 

-4.04352E-19 


31 

0- 

-1,617826-18 


ORIENTATIONS 

STRESS 

ORIENT ATIONS 

TRUSS 

FIBERS IN 

TRUSS 

FIBERS IN 


THU5S TRUSS 


MAXIMUM 

Stress 

*•000006+04 
4.00000E+04 
4*000006+04 
4*000006+04 
4. 000 Q 06 + 04 
4.000006+04 
4.0OOQOE+O4 
4,000006+04 
4, 000006+04 
4* 000006 + 04 
4,000006+04 
4-000006+04 
4# 000006*04 
4 , 000006*04 
4.00000E+04 
4.00000E+04 
*• O0OOOE+O4 
4* O00OOE+ 04 
4* 000006+ 04 
4- 000006+04 
4,000006+04 
4* 000006+ 04 
4, 0OOQOE + O4 
4* OOOOOE+04 
4. OOOOOE+04 
4iOOOOOE*04 
4. OOOOOE + 04 
4*000006+04 
4i 000006+04 
4, OOOOOE+04 
4iOOOOOE+04 
4, OOOOOE+04 
4, 000006+04 
4.000006+04 
4. 000006 + 04 
4-000Q0E+04 
4.00000E*04 
4.000006+04 
4.000QQE+04 
4.00000E+04 


STRFSS 



REPRODUCIBILITY of rnv 
ORIGINAL PAGE IS 


1 

1 

1 .3674 Jt+03 

1 

2 

1.297356*03 

1 

3 

) . 1 56336* 03 

1 

4 

9.435706*02 

1 

5 

6 * 591 1 St ♦ 02 

1 

6 

3.111436*02 

1 

7 

-6 .469036*01 

1 

8 

-3.40819E* Q2 

2 

I 

1 # 367736* 03 

2 

2 

1 * 300006*03 

2 

3 

] *16395t*03 

2 

4 

9.585526*02 

2 

5 

ft. 63953fc*02 

? 

6 

3.480336*02 

2 

7 

■1 *477836* 0 1 

? 

S 

-2.890696*02 

3 

l 

1.367736*03 

3 

2 

1.3000BE*O3 

3 

3 

1.163956*03 

3 

4 

9.585526*02 

3 

5 

6.839536*02 

3 

6 

3. 48033E* 02 

3 

7 

-1.477826*01 

3 

e 

-2.890696*02 

4 

1 

1 , 3680-16*03 

4 

2 

1.302816*03 

4 

3 

1*1 7 156E* 03 

4 

4 

9.73533E*02 

4 

5 

7.087686*02 

4 

6 

3.849236*02 

4 

7 

3.513366*01 

4 

8 

-2.293186*02 

5 

1 

-2.701 18E* 02 

5 

2 

-2.749746*02 

5 

3 

-2.847476*02 

S 

4 

-2.99490E*02 

5 

5 

-3.191986*02 

5 

6 

-3.433046*02 

5 

7 

-3.693356*02 

5 

6 

-3.890136*02 

SOLUTION FOR STRESS 

converges IN 

2 ITERATIONS 






direction 

APPLIED LOAD 


STRAIN 



11 

0* 


-1.49926E-04 



22 

0* 


-1.0306BE-04 



33 

1.100006*01 


5.25434E-04 



12 

0. 


9.369656-13 



23 

0* 


-4*409996-18 



31 

o. 


-1.460456-17 



ORIENTATIONS 

STRESS 


ORIENTATIONS 

STRESS 

TRUSS 

FIBERS IN 

TRUSS 


fibers IN 



TRUSS 



TRUSS 


1 

1 

1 .496116*04 

1 

2 

1.419656*04 

1 

3 

1.265706*04 

1 

4 

1.033236*04 

1 

5 

7.220486*03 

1 

6 

3.407046*03 

1 

7 

-7.172436*02 

1 

a 

-3.63735E*03 

2 

1 

1.496446*0* 

2 

2 

1.422636*04 

2 

3 

1.274036*04 

? 

4 

1.049636*04 

2 

5 

7.492026*03 

2 

6 

3*81 307E* 03 

2 

7 

-J .681086*02 

2 

8 

-3.180106*03 

3 

1 

1 ,496446*04 

3 

2 

1.422636*04 

3 

3 

1.274036*04 

3 

4 

1.049636*04 

3 

5 

7.492826*03 

3 

6 

3.813076*03 

3 

7 

-1.681006*02 

3 

8 

-3.18010E*03 

4 

1 

1 .49677E* 04 

4 

2 

1.425626*04 

4 

3 

1 .202366*04 

4 

4 

1.066046*04 

4 

5 

7 .765136*03 

4 

6 

4.218216*03 

4 

7 

3.81 007t*02 

4 

8 

-2.522816*03 

5 

l 

-2.970796*03 

5 

2 

-3.024216*03 

5 

3 

-3.131716*03 

5 

4 

-3.293896*03 

5 

5 

-3.510686*03 

5 

6 

-3.775846*03 

5 

7 

-4.06217fc*03 

5 

8 

-4*278636*03 


SOtOtTOw FOR STRESS CONVERSES IN ? ITERATIONS 



direction 


APPLIED LOAD 


STRAIN 



11 

o. 


-2.86225E-04 



22 

0. 


-l^GTSOE-O* 



33 

2* 10000E+01 


1.00109E-03 



12 

o. 


1.T6031EM2 



23 

0* 


-T,83209E-1B 



31 

0. 


-2.81076E-17 



ORIENTATIONS 

STRESS 


ORIENTATIONS 

STRESS 

TRUSS 

FIBERS IN 

TRUSS 


FIBERS IN 



TRUSS 



TRUSS 


1 

1 

?• B4103E+04 

1 

2 

2»69622E*04 

1 

3 

?,404S?E*04 

1 

4 

1.96365E+0* 

1 

5 

W372B1E*04 

1 

6 

6.47T92E*03 

1 

7 

•l*3B003E* 03 

1 

8 

«.?.3Z847E*03 

2 

1 

2.B416GE*04 

2 

2 

Z*7018BE*<H 

2 

3 

2.42033E*04 

2 

4 

1 • 994B3E*Q4 

2 

5 

1.42464E404 

2 

fi 

7 »25045E*Q3 

2 

7 

-3*31490E*02 

2 

8 

-6.07173E*03 

3 

1 

2.84166E+04 

3 

2 

2,70188E*0A 

3 

3 

2*42033E*OA 

3 

4 

1-99483E*0A 

3 

5 

1.424S4E*0A 

3 

b 

7,25045£*03 

3 

7 

-3.31490E*02 

3 

8 

-6.07173E+03 

4 

1 

2.B4228fc*04 

4 

2 

2.70755E*04 

4 

3 

2»43614t*04 

4 

4 

2.02601E+0A 

4 

5 

U47646E*04 

4 

b 

0» 0226BE*03 

4 

7 

7* 16983E *02 

4 

8 

-4.81687E*03 

5 

1 

-5.67053E*03 

S 

2 

-S*77ZS2E*03 

5 

3 

•5»97775£+03 

5 

4 

-6.207376*03 

5 

5 

-6.7012SE+03 

5 

b 

-7.20T45E*03 

5 

T 

*T»T540TE* 03 

5 

8 

**B. 1G728E*03 


SOLUTION POP STRESS CONVERGES IN 2 ITERATIONS 



direction 

APPLIED LOAD 


STRAIN 



11 

0* 


-4.22527E-04 



22 

0* 


-2.9O417E-04 



33 

3.10000E*01 


1.4T404E-O3 



12 

0- 


2.63917E-12 



23 

0. 


-1.Z881RE-17 



31 

0. 


-5.00643E-17 



orientations 

STRESS 


orientations 

STRESS 

TRUSS 

FIBERS IN 

TRUSS 


FIHERS IN 



TRUSS 



TRUSS 


l 

1 

4 • 17181E*04 

1 

2 

3.95972E*04 

1 

3 

3.53228E404 

1 

4 

2.B85T4E*04 

1 

5 

2»0lB2e»E*04 

1 

6 

9.52159E*03 

1 

7 

-2.05301E*03 

1 

8 

-1.08182E*04 

2 

1 

4.17273E*04 

2 

2 

3 ff 96BO3F*04 

2 

3 

3.555SOt*04 

2 

4 

2.93150E*O4 

2 

5 

2.09454E+O4 

2 

b 

1 • 06604E*04 

2 

7 

-5.0487GE402 

2 

8 

-8.9*395E*03 

3 

1 

4*1 7273E* 04 

3 

2 

3*9ftB03E*04 

3 

3 

3.5S550E*04 

3 

4 

2*93158E*04 

3 

5 

2*09454t*04 

3 

6 

1 • 06S04E* 04 

3 

7 

-5.04B76t*02 

3 

B 

-B.96395E*Q3 

4 


4.l7364t*0A 

4 

2 

3.97G34F*04 

4 

3 

3.S7872E *04 

4 

4 

2.97740E*04 

4 

S 

?_ } Titans «ni 

A. 

a 

i 1 Ttatt t. c * tlA 



4 

5 
5 
5 
5 

MEMBER FAILURE! 
MEMBER FAILURE! 
MEMBER FAILURE! 
MEMBER FAILURE! 


7 

1 

3 

S 

7 

TRUSS 1 FIBER 

TRUSS 2 FIBER 

TRUSS 3 FIBER 

TRUSS 4 FIBER 


1.04311E+03 

-R*3693JE+D3 

•B.8228fe£*03 

-9.09O88E>O3 

-1«14450E*04 

1 

I 

1 


4 

5 
A 
A 
A 


ft 

2 

4 

6 

6 


11 14flE*03 
-A.519B9E*03 
•9.27992E*03 
-1*0A3B1E*04 
-p!.20549E*04 


SOLUTION FOR STRESS CONVERGES IN 2 ITERATIONS 



i*r- 

c : 



direction 




APPLIED LOAD 


STRAIN 


1 1 




0* 


-5,5B636E-04 


22 




0* 


-3.B393SE-04 


33 




2 t96957E*0 1 


l*95566E-03 


12 




0. 


3 t 4B834£»l 2 


23 




0. 


•1 #2G954E*l7 


31 




0« 


•4*246S7E^1 7 


ORIENTATIONS 




STRESS 


orientations 

TRUSS 

FIBERS IN 




TRUSS 


FIBERS IN 


TRUSS 






TRUSS 

1 

1 




0* 

1 

2 

1 

3 




4* 67136E* 04 

1 

4 

1 

5 




Z.67383E*04 

1 

6 

1 

7 




-2.68378E+03 

1 

8 

2 

1 




0. 

2 

2 

2 

3 




4,70190E*04 

2 

4 

2 

5 




2.7743Bfc*04 

? 

6 

2 

7 




•6.36594t*02 

? 

B 

3 

l 




0* 

3 

2 

3 

3 




4,70 190E*OA 

3 

4 

3 

5 




2,77438E*04 

3 

6 

3 

7 




-fc*36594E*02 

3 

e 

4 

1 




0. 

4 

2 

4 

3 




4*73243E*0A 

4 

4 

4 

5 




2,fl74BBt*04 

4 

6 

4 

7 




1 * 41031E*03 

4 

6 

5 

1 




•1*1 0633E ♦ 04 

5 

2 

5 

3 




*1 . 16630E+0 A 

5 

4 

5 

5 




-1*30750E*04 

5 

ft 

5 

7 




-1*51297E404 

5 

0 

MEMBER 

FAILURE! TRUSS 

1 

FIBER 

2 




MEMBER 

FAILURE! TRUSS 

2 

FIBER 

2 




MEMBER 

FAILURE! TRUSS 

3 

FIBER 

2 




MEMBER 

FAILURE! TRUSS 

4 

FIBER 

2 





STRESS 


S*2344?E*04 
3*8l 89?E 4 04 
1 • 26A44E*04 

-1,42946E*04 
5.24539E*04 
3* BT925E*04 
1*41479E*04 
-1*1845BE*04 

5*24539E*0* 
3 • B7925E^04 
1.41479E*04 
•1.18458E404 
S.2563oE*Q4 
3.93957E*04 
1, 56504*404 
•9*39659E»03 
-1.12624E*04 
-1.22673E*04 
“1 * 40829E*04 
-1.59360E-M)4 


SOLUTION FOR STRESS CONVERGES IN 3 ITERATIONS 


DIRECTION 

APPLIED load 

strain 

11 

0# 

•0 # 22897E — 04 

22 

o. 

-5.BS390E-04 

33 

2»V69S7E*f>l 

2*99701 F*03 

12 

0* 

S.13ft54E-l2 


ft . 

* 1 ^TAQQtk* Tf 



31 

ORIENTATIONS 

TRUSS FIBERS IN 

TRUSS 


0 - 

5TRES5 


1 1 

1 3 

1 5 

1 7 

2 1 

2 3 

2 5 

2 7 

3 1 

3 3 

3 5 

3 7 

* 1 

4 3 

4 5 

4 7 

5 1 

5 3 

5 5 

S 7 

MEMBER FAILURE I TRUSS 


0 * 

7 • 127526*04 
4 » 12136E+-04 
~3.33460E*03 

0 . 

7,1?200E*04 
4,2685 IE*04 
•3, 171 30fc>02 

0 , 

7.17200E*0* 

4,Z6851t#04 

*►3, 171 30E*02 

0. 

7,21 646E + 04 
4, 41555E* 04 
2,6989U*03 
•1 , 62891E+ 04 
•1 , 71727E+04 
•1 ,92531E*C4 
-2, 22801 t+04 

4 FIBER 3 


solution FOR STRESS CONVERGES IN 3 iterations 


direction 


applied load 


0 * 

0 * 

2.9695TE*0J 

0 ' 

o. 

0 . 

STRESS 


1 I 

1 3 

) 5 

1 7 

2 1 

2 3 

2 5 

2 7 

3 1 

3 3 

3 5 

3 7 

4 1 

* 4 3 

cn 4 5 

4 7 

5 1 

5 3 

5 5 

5 7 


MFMRFR PiTI URF i 


0* 

7.97I29E+04 
4.59116L+04 
"4 ,2B974fc*o3 
0 . 

R# 03484E* 04 
4* 80 176t* 04 
3* 92926E* 01 

0 » 

8* 03484E + 04 
4,801 766>04 
3,92928E*01 

0. 

0, 

5, 0121 4E 404 
4 ,36423E*03 
-1.66427E404 
-1 ,791 02t 404 
-2,08947fc>04 
- 2 . 52368 E 404 


TRII^S ? 


PTRFH 1 


I #7 < 4.1 

-8. 12669E-17 



ORIENTATIONS 

STRESS 

TRUSS 

FIBERS in 



TRUSS 


1 

2 

0. 

1 

4 

5* 84641 E* 04 

1 

6 

1 • 991 50£* 04 

1 

8 

-Zt09701E*04 

2 

2 

0. 

2 

4 

5t93446E*04 

? 

6 

2* 21 229E 4 04 

2 

8 

• l * 73622E* 04 

3 

2 

0, 

3 

4 

5*93446E*04 

3 

6 

2.21ZZ9E*04 

3 

8 

«• 1 * 7362ZE* 04 

4 

2 

0, 

4 

4 

6, O2240E*O4 

4 

6 

2»43204E* 04 

4 

8 

«*l , 37533E* 04 

5 

2 

*1 ,65825E*04 

5 

4 

-1,80630E*04 

5 

6 

-2,07086E404 

5 

8 

*2*34681 E* 04 


STRAIN 

W 9,465B4E*0* 

-5,77144E-04 

3,36218E*03 

5.02068E-12 

w|,37373E-l7 

-7,05499E-]7 



ORIENTATIONS 

STRESS 

TRUSS 

FIBERS IN 



TRUSS 


1 

2 

0* 

1 

4 

8»53154E*Q4 

1 

6 

2*19268E*04 

1 

8 

-2*41791E*04 

2 

2 

0* 

2 

4 

6,65744E*04 

? 

6 

2*50905E* 04 

2 

8 

-1*90039E*04 

3 

Z 

0, 

3 

4 

6* 65744E+04 

3 

6 

2«50905E*04 

3 

8 

-1 .900395*04 

4 

2 

0* 

4 

4 

6.78326E»04 

4 

6 

2 ,8?493E*04 

4 

0 

• 1 • 38266E*04 

5 

2 

-1*70635E*04 

5 

4 

-1.91875E-MJ4 

5 

b 

• 2* 29828F* 04 

5 

B 

»?* 69407E ♦ 04 



SPACE PROGRAM LISTING 


APPENDIX D 

The following pages present a listing of the CDC 6600 
version of the SPACE program. The program is written for 
a "RUN" compiler, and requires 70K of core to compile and 
execute. The compilation of the program and execution of 
the sample case given in Appendix C required less than 13 
system seconds of time. 
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HUN VERSION 2*3 — PSR LEVEL 363' 


08/?2/74 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

000003 

000003 

000003 

1>0I>003 


000003 

000003 

000003 

000003 

000003 

000003 

000003 

000003 

000003 

c 

c 

000003 

000012 

OOOOU 

0000 ?* 

000033 

000060 

000101 

000114 

C 

C 

c 

000123 
0001 ?5 
uom ?*j 

ooni?7 

000130 

000132 

00013* 


PROGRAM MAIN ( INPUT . OUTPUT *TAPE5* INPUT *TAPE6«0UTPU I] 


4 SPACE 

« 

* VERS X MOD 0 


COMMON /CINPUT/ADtBD*00*NFl f L*EY*eCtA*VMQD*VTl*V*5GM 
C0MM0N/L0C2/5T0 
DIMENSION TEMPO) 


DIMENSION 

1 

2 

3 

4 

5 


DIMENSION 

1 

DIMENSION 

namelist/data/ 

1 

DATA 


EY ( 5* 30) * A ( 5*30 ) *EC ( 5*30>*BT< 5*30>*NB< 5*30)* 
U 5*30)»ALr 5) *SOO < 6 * 1 > ¥TRIG U *6) tTRMO < 6.1>* 

EV ( 6« 1 ) * C ( ?1* 5*3q) *OT6t6) #DIRFU*1> *SGMl 5*30)* 
ETA ( 5*30>*EP< 5*30>*RAT< 5*30>*5GF< 5*30>* 

SGP ( 5*30) iRI 6* 5*30) *E(6*6) *Z<6*6) * DE ( 5*30)* 
DM ( 5*30) * 5TOH5 ( 5*30)*$M( 6)*VT< 5*iO> 

5LN£6*1) *S(_Pt&*l) *SGC( 5*30)*VTK 5*30)* 

FRAC ( 5t30)*TMPFL( 5*30*40) 

II (40) * JJ(4(>) 

NIT*SEP5»SUPB*NOR*NTR*AD»0D*DD*$QM*RAT*EC*EY*Li 
VTI*TEMP t ALl > HA*NB*AL*SGOiNlNCtMXCALtSMLT f RT 
MXCAL*NINC*NlT*SMLT*SEP$*SUPa 

Jnn.R.i 


.1 rtC * nc J 


REAL N0,L 

INTEGER STO (5*30) 

FQUIVALENCE (OE*EP> 

data Pl/3,1415q?6EO/*EXPAHG/l73.0E0/ 

LOGICAL PSBT*MS1NGD*IFAIL( 5*30) 


DSECCX) ■ 1.0E0/ COS<X) 

RAD(X) « K/57-29570 

DCOTHU) « *.0E 00/ TANH(X) 
um * SORT C X) *P| /2« 

COEL ( C • A ) * C*A«M1.-IP1/2*)*<1.- CO$( U(A)))/< U(A)* 

1 SIN ( U(A))))/n # -A) 

TDEL 1 C* A) ■ C # A* 1 1 *- IPI/2*) / l UCA)*OCOTH< U(A>/2,)))/ 
1 (]*♦*) 

EVAI<Y*A*C*D) * A*Y/(A*4t* (2# # C*D*D**2) ) 

F ( X) ■ 1.00E 00 - EXP(-X) 


INPUT 

AD*1 • D 

RD= i . 0 

DO*! < 0 
I C5 c C 

1 TCS»IC5*1 
READ ( 5* DaT A) 

IF < EOF *5) 33*32 


§ 

£ § 




RUN VEflSTON 2.3 — PSR LEVEL 363' 


MAIN 


OB/22/74 


000137 


□00141 

000143 

000144 


00015? 

000154 

000155 

000160 

000164 

000165 

000166 

000167 

000170 

000171 

000172 

000174 

000177 

000200 

00020? 

U00203 


000205 

000206 

000207 

000211 

000227 

000242 

000242 

000251 

000266 

000272 

000274 


000300 


C INPUT PER CASE 

C 

c CONTROL PARAMETERS 

32 NORP* NOR * 1 

c 

C RESET MATRIX STO TO ZERO 

C 

00 19 KDF* 1*5 
00 19 KDZ* 1,30 
19 STO < KDF » KDZ ) *0 

C MODEL GEOMETRY 

C STRESS LIMIT 

C FIBER AND TRUSS PROPERTIES 

C ANGULAR ORIENTATIONS 

C INTIAL LOADING 

C 
C 

C INITIAL COMPUTATIONS AND ASSIGNMENTS 

DO 440 I.ltNORP 
DO 440 sJ*l * NTR 
JFAIL 1 1 * J) = .FALSE. 

440 CONTINUE 

MSINGD * .FALSE. 

INC ■ 0 

IT * 0 

TTOT * n 

MEMB - 1 

VMOD ■ O.OE 00 

KNT « NORP*NTR 

V ■ AD*RO*DD 

DO 480 1 « 1 t 6 

SM(I) * SMLT*SGO£I*l> 

SLN ( 1 * 1 ) = 0 . 0 E 00 
480 CONTINUE 

c 

c 

DO 910 I a 1 * NURP 

DO 910 Jsl *NTR 

IF t INP.EQ. 2 ) GO TO 090 

> UI •J)*DSEC1RAD tRTIIfjm/ t2,0*NB(I .JI*HAT CItJ) > 
ETA 1 1 * J) ■ 4**NB £ I . J) **2*DD/ ( L < l * J) * T AN ( RAO < BT <1 * J) ))) 

090 CONTINUE 

A<I,U> * Pl*DM(ItJ>#* 2 / 4 . 0 EO 

VT(t»J) « 2 .*NR(liJ)*A<I*J>/( SINCRAUiBTUf J)))»AO«BD) 

FRACUtJ) * VTI UtJ)/VT<IiJ) 

VMOD * VMOU ♦ VT (I, J) 

910 CONTINUE 
C 

C PRINT OUT INPUT INFORMATION 

CALL INPUT <NUR*NTR*ICS» AL » RT ) 

C 

C 

C OEGREE TO RADIAN CONVERSION 

DO 912 I» 1 * NOR 
AL(I> * RAOiALUn 


000304 

000306 
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000312 
000315 
00031 6 
000317 
000324 


000331 

000331 

000333 

000334 

000337 

0003*7 


000354 

000354 

00035ft 

000357 

000352 

000355 


000371 

000372 

000375 

000376 

000377 

000400 

000400 

000402 

000403 

000406 

000413 

000414 

000414 

000416 

000417 

000422 
000433 
000436 
000443 
000444 
000445 
000447 
000450 
000454 
000455 
000463 
000464 
^ 000467 


912 CONTINUE 

DO 913 1*1 *NORP 
00 913 J*1*NTR 
RT<J*J> ■ RACHBT(IiJ) > 

913 CONTINUE 
C 

c RETURN TO 915 TO INCREMENT LOADING 
915 CONTINUE 
C E* AT ALPHA * 0 

00 917 1*1* NORP 
no 917 J*liNTR 
IF ( IFAIL ( I * JH GO TO 917 

OEUfU) ■ Er<ItJ>/(8.0Eo* EC<l*J)** 2 * 1 . 0 E 0 ) 

917 CONTINUE 

c RETURN TO 920 FOR NEXT ITERATIVE STEP 

C INITIALIZE MATRICES 0 ♦ E TO ZERO 

920 CONTINUE 

00 950 M*J *6 
DO 950 N*lf6 
D<M.N> * 0 « OE 00 

FtMiNl * O.OEO 
950 CONTINUE 
C 
C 

DO 1005 1*1*6 
TRMOtlfl) * V*S60Uil> 

1005 CONTINUE 

C EVALUATE MATRICES D ♦ E* FIRST FOR +BETa ♦ SECONDLY FUR - 0ETA 

PSBT « .TRUE. 

00 TO 1020 
1010 CONTINUE 

00 1015 I*l*NORP 
00 1015 J* 1 *NTR 
RT(I*J) ■ -BT(I*J) 

1015 CONTINUE 

PSBT ■ .FALSE, 

1020 CONTINUE 
C 

no 1750 1*1* NOR 
DO 1750 J* 1 *NTR 
IF ( IFAIL (I * J> 1 GO TO 1750 

C SET PRELIMINARY VALUE OF MATRIX ELEMENTS 

PP * ETA(I* J)*nM(I*J)*RAT<I* J)»AlI*JT/2.*FRACtI*J) 

SNA * SIN(ALII) ) 

SN2A * SIN<2t*AL(I) ) 

SNA? ® SNA**2 

SNA3 * $NA2*SNA 

SN2A2 * SN2A#*2 

SNA 4 * 5NA2##2 

SNR2 * SIN <BT ( I * J) )**2 

SNB4 ■ SNS2**2 

5N2B * SINT 2.*BT ( 1 1 J) ) 

SN2B2 * SN28*«2 

CSA * COS<AL<I>) 

CSA? * CSA**2 


g? Q 

,e & 

^ Cy 

se 

o 

k 



LP 

O 
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U0ft4?0 

000471 

000472 

000476 


000477 

000506 

000515 

000525 

000535 

000547 

000557 

000567 

000575 

000605 

000615 

000627 

000635 

000644 

000654 

000665 

000674 

000704 

000715 

000724 

000733 

000742 


000747 

000750 

000751 

000752 

000754 

000755 

000756 

000761 

000774 


001 005 
001007 

ooion 

001012 

001017 

001023 
001025 
001 026 
001031 
001042 
001046 
001047 
001055 


CSA3 ■ CSA2*CSA 

CSA4 * CSA2**2 

CSB2 * COS (BT £ I * J) ) **2 

CSB4 « CSB2**2 


ELEMENTS OF 

D MATRIX 

Cl 

1 $ I * J> 

a 

PP*CSA4*C504 

Cl 

2« I f J> 

a 

PP* CSA2*SNA2*CSB4 

Cl 

3 a I * J > 

a 

PP*CSA2*CSB2*SnB2 

Cl 

4 * I * J) 

■ 

-PP*CSA2*SN2A*CSB4 

Cl 

5 * I » J) 

a 

-PP*C5A2*SNA*CS92*SN2B 

c< 

6* I * J) 

a 

-PP* CSA3*CSB2*SN26 

Cl 

7»lt J) 

m 

PP*CS04*SNA4 

Cl 


a 

PP*CS02*5NB2*SNAZ 

Cl 

9t I « J) 

a 

-PP*C5B4*5NAZ*SNZA 

Cl 

lOfliJ) 

a 

-PP«5NA3*CS02*SN2B 

c< 

11*1 * J) 

a 

-PP*CSA*SNA2*CSB2*SN2B 

cl 

12*1* J) 

a 

PP*5NB4 

Cl 

1 3f I « J) 

a 

-PP*SN2A*CSB2*SNB2 

c< 

1 4f I ♦ J) 

a 

-PP*SNA*5NRZ*5N2B 

Cl 

15* I* J> 

a 

-PP*CSA*SNR2*5N2B 

Cl 

16* 1 * J> 

a 

PP*SN2AZ*C5B4 

cl 

17* I * Jl 

a 

PP# SNA*SN2A*CSR2*SNZB 

Cl 

1 8 * I * J) 

a 

PP*CSA*SN2A*CSB2*SNZfl 

Cl 

1 9 1 1 * J > 

a 

PP*SNA2*SN2B2 

Cl 

?0i 1* J) 

a 

PP*5NA*C5A*SN?R2 

cl 

21*1 * J) 

a 

PP*CSA2*5N2B2 


1750 CONTINUE 
C 

c FORM UPPEH TRIANGLE OF D MATRIX 

MC * 0 

00 1600 Mai >6 

00 1600 N*M » 6 

MC * MC ♦ 1 

no 1600 I ■ 1 • NOR 

no 1800 Jal »NTR 

IE l IFAIL f T * J) ) 00 TO 1800 

0 <M*N) « DlM.N) * C( MC f I f J) *OE 1 1 » J) 

\800 CONTINUE 
C 

c FORM SYMMETRIC ELEMENTS OF 0 
DO 1850 M* 2'6 
LO a M*t 
DO 1050 N* 1 *LD 
0 (M,N) « U(N,M) 

1850 CONTINUE 

c 

1 a NORP 

no 1900 JaltNTR 
IFdFAlKltJH GO TO 1900 

PP a ETA<I»J)#OM(I*J)*RATn »J)*AU *J) / 2 »*FRAC(It J) 

SNB 2 * SIN(BT ( I. J> )•*? 

SNB 4 a 5 NB 2**2 

SN 2 R a SIN( 2 .* 8 T(I*J)) 

SN 2 R 2 » SN 2 B **2 
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001056 

0010*2 


001063 

001071 

001076 

001102 

001106 

001113 

001117 


00112 ? 
001 1 23 
001126 
001134 
001143 
001151 
001160 
001166 
001175 


001200 

001201 

001203 


001204 

001206 

001207 

001210 

001213 


001220 

001220 

001222 

001223 

00123? 

001236 


001237 


001246 

001250 

001251 

001260 

001261 

001261 

001277 


CSB2 * C0S<BTCI«J)>**2 

CSP4 * CSB2**2 

C 

C ELEMENTS OF E MATRIX 

G( lilij) * PP* CSB4 
G ( 2*1 »J) « PP*C5B2*SNB2 

Gl 3 ♦ I * J) » -PP*CSB2*SN2B 
G< 4*1 • J) » PP*SNB4 
G ( 5* I * J J * -PP*SNB2*SN20 

G< 6.1*J) * PP*SN2B2 
1900 CONTINUE 


c 

c 

FORM UPPER TRIANGLE 
DO 1910 J*l*NTR 

OF E MATRIX 


IFllFAlL<1 

:*j> ) gc 

1 TO 

1910 


E<1*11 • 

E ( 1 * 1 1 

♦ G ( 

1 *I*J)*DE(I* JJ 


E ( 1 9 2 ) * 

EM *2) 

* Gf 

2 * I * J) *DE ( I • J ) 


?U»T) « 

E ( 1 * 41 

♦ G( 

3 * I * J) *UE ( 1 1 J) 


E (2*2) * 

E ( 2 * 2 ) 

♦ G ( 

4* I * J) *DE < I * J> 


F (2*4) « 

E ( 2 9 4 ) 

♦ Gt 

5* I * J) # DE ( 1 1 J) 

1910 

E (4 * 4) * 
CONTINUE 

E ( 4 *4 ) 

♦ G( 

6*1*J)«0EU*J> 

c 

c 

FORM SYMMETRIC 

ELEMENTS OF E 


E (2f 1 ) ■ EUt?) 

E(4*l) * E 1 1 1 * > 

EU*2> • EIZ.4) 

C 

C CHECK FOP *BETA OR -BETA 

TF (PSBT 1 GO TO 1010 
00 1915 I* 1 f NORP 
00 1915 JBltNTR 
RTU.J) ■ -RT(1,J1 
1915 CONTINUE 
C 

c form sum of matrices d ♦ e 

1917 CONTINUE 

DO 1920 M«li6 
00 1920 N*1 1 6 
Z(M*N> * 0 <M * N) * E(M*N) 

1920 CONTINUE 

lF(MSlNGP) GO TO 1967 
C 
C 

C INVERT RESULTANT MATRIX 

CALL INVRTD(Zf 6 »6tOET*l.E»15* IRANKtl. 00F-50) 
C 

c CHECK FOR SINGULAR MATRIX 

IF { IRANK,E0.6> GO TO 1968 
MSINGO ■ .TRUE. 

WRl TE <6« 97 1 0 > JRANK,OET 
GO TO 1917 
1967 CONTINUE 

WRITE <6*97 15) ( <Z < K 1 .LI) *11*1 »6> .Kl = l *6) 
MSINGO * .FALSE* 
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001300 

001301 


001301 

001301 

001303 

001304 

001307 

001311 

001313 

001314 

001321 

001324 

001325 

001325 

001334 

001340 

001345 


001347 

001351 

001353 

001354 

001356 

001360 

001361 


001361 

001362 

001367 

001367 

001400 

001402 

001407 

001410 

001421 

001423 

001424 

001425 

001435 

00144? 

001442 

001445 

001460 

001470 

001476 

001500 

001501 

0015^1 


1968 

C 

C 

2000 


2050 


C 

C 


2080 

C 

c 

c 


3110 

C 


3115 


3201 

C 


3202 

3525 


60 TO 7000 
CONTINUE 


STRESS CALCULATIONS 
CONTINUE 

DO 4000 1*1* NORP 
DO 4000 J«l*NTR 
IF<IfAILn*Jn GO TO 4000 
IF(I*EQ.NORP> GO TO 2050 
SNA * SlNtALUM 

5NA2 ■ SNA**? 

SN2A * SIN(2,#ALU>> 

eSA • cost AH I)) 

CSA2 * CSA##2 
CONTINUE 

5N2B ■ SIN(2*«8TU*jn 

SN02 * SIN ( 0T < 1 1 J ) ) **? 

CS&2 ■ COStBTUr J>>**2 

IF (I.EQ.NORP1 GO TO 2000 


stress-strain factors 
TRIG ( 1 i 1 ) * CSA2ACSB2 
TRIG ( 1 *21 « CSB2*SNA2 
TRlGtl t 3) ■ 5NB2 
TRIG<1*4J ■ -SN2A*CSB2 
TRIG(1*51 » -SNA*SN?9 
TRIG 1 1 *61 ■ -C5A*SN?B 
CONTINUE 


STORE STRESS VALUE FROM PREVIOUS ITERATION 
IF(IT.EQ.O) GO TO 3110 
SGPUtJ) • SGF(1*J> 

CONTINUE 

CALCULATE STRAIN values 

CALL MXMULHZ* TRM0*EV*6*6* 1 *6*6*1 > 

DO 3115 ITMP*1,3 

EV(ITMP»11 * EV(ITMP*11 * *5*ALPHA*TEMP £ I TMP) 

CONTINUE 

CALL HXMULDtTRIG*EV,DIRF, 1*6, It 1.6*1) 

COMO ■ DIRF U • 1 ) 

TRIG ( 1 * 5> ■ -TRIG ( 1 *5) 

TRIG ( 1 * 6 > « -TRIG <1*61 

CALL MXMULO { TRIG*EV iDlRf»lv6«l*]t6»l) 

IF t ABS < COMD) *LE ■ ABS(DIRFMflM) GO TO 3201 
D I RF Civil * COMD 
IF (I.LT*nO«P) GO TO 3525 
CALCULATE STRESS VALUES 

SGF 1 1 * J) ■ EP(I*J)*(EV(l*n*CSR2*EVtZ*ll*SNR2-EV(4.1»ASN2B) 
COMD * EPU * J)*(EV< 1 *11*050? ♦£ V ( 2» 1 ) *SNR2 *E V 14* 1 J *SN?R ) 

TF C ABS < COMD) »LE ■ ABStSGF 5 1 * J) ) > GO 10 3202 
SGF 1 1 * J) * COMD 
GO TO 3550 

continue 

SGFtlfJ) ■ EP t!*U) ADIRF (1,1 ) 
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00150? 
00 1 507 
001522 
001523 

001543 

001563 


001570 

001571 

00157? 

0015^3 

001573 

001575 

001576 

001601 

001606 

001617 

001621 

001621 

001624 

001627 

001630 

00163? 

001633 

001634 

001635 

001636 

001643 

001644 


001653 

001653 

001657 


001660 

001660 

001664 


001665 

001665 

001673 
u, 00 1676 


3550 CONTINUE 

ALP * l6. # SGF(IfJ)/<Pl**2*PV U,J> ) • (RAT < I* J>> **2 
A8SA « ABS(ALP) 

IFISGF<1, J) .LE.O.OE 003 EP<I*J> » EV AL ( EY ( I * J> * ABSA • EC l 1 * J> • 
1 CUtL < EC ( I * J) l ABSA > > 

IP(SGFU»J1.GT.0.0E 001 EPU»J> * EV AL < EY < I ♦ J) » ABSA • EC t I • J) t 
1 TUEL lEC < I * J) t ABSA) ) 

4000 CONTINUE 
C 

C STRESS CONVERGENCE TEST 
IF(IT.NE.O) GO TO 4010 
IT - 1 
GO TO 92o 
4010 CONTINUE 

00 5000 1*1 *NORP 
00 5000 J«1.NTR 
IF ( IFA1L C I * J 3 ) GO TO 5000 

DSG - ABStSGFU, Jl-SGPCI* Jl) 

C CHECK FOR ZERO STRESS OR RELATIVE STRESS ,LT. 1.0-06 

IFMSGP II* J) .EO.0-3 -OR. <O5G.LT.l.E-06M GO TO 4B0Q 
OSG ■ ARS ( DSG/SGP ( I * J) ) 

4RQ0 CONTINUE 

IF (QSG.lt , 5EPS) GO TO 5000 
IFtOSG.GT. 5UPB) GO TO 4900 
C STRESS ITERATION 

IT * IT * 1 
IF(IT.LT.NIT) go TO 920 

C NO CONVERGENCES WITHIN ITERATION LIMIT 

ICON » 1 
GO TO 5100 
C DIVERGENCE 
4900 ICON ■ ? 

GO TO 5100 
5000 CONTINUE 
C NORMAL CONVERGENCE 
ICON ■ 3 

5100 GO TO ( 5?00t5250f 5300 ) * ICON 
C 

c NON-CONVERGENCF DUMP 

5?00 CONTINUE 

CALL COUMP(NURP*NTRi IC0N*1T*SGF,SGPJ 
GO TO 7 000 
C 

C DIVERGENCE oump 
5250 CONTINUE 

CALL COUMP (NORP *NTR* 1 CON , I T ,SGF *5GP ) 

GO TO 7000 
C 

C CONVERGED SOLUTION ANALYSIS 

C 

5300 CONTINUE 

WRITE (6*97201 IT 

C CHECK * INCRIMENTATION LIMIT 

IF t INC.GE.NINC> GO TO 7000 
INC * INC ♦ 1 








& 



U1 

4* 
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001677 


001700 

001702 

001703 

001705 

001706 


001710 

001710 

001712 

001713 

001716 

001730 

001731 

001762 


001776 


002026 

002033 

002033 


002033 

002037 

002061 


002063 

002063 

002066 

002066 

002067 

002052 

002061 

002062 

002071 

002072 

002076 

002101 

002102 

002116 


IT ■ 0 

C 

C STORE LOADINGS FOR PRESENT ANO PREVIOUS lNCftlMENT AT IONS 

|FtMEH9.N€*l) GO TO 5955 
DO 5950 1-1*6 

SLP(Ifl) » 5LN < I « 1 > 

SLN ( 1*1) ■ SGO ( I * 1 > 

5950 CONTINUE 


5955 


COMPUTE MAXIMUM STRESS COMPARISON TERM FOR EACH MEMBER 
CONTINUE 

00 5960 T»1*N0RP 

00 5960 J»1»NTR 

IF < 1FA1L ( I » J) ) GO TO 5960 

ALP - 16**SGF (I«J)/CPI**2»EY(X'J) ) * ( RAT t It J> ) **2 

ALB ■ ABS(ALP) 

TRM 2 * PI*ET (I* J)*eC<I«J)*U(ALB)/RATII»J)«*z 
lF(5GF<ItJ).LE.0*0E 00) SGC 1 1 * J) « SGF(ItJ) * 


IF ( SGF < 1 1 J) *GT*Q*OE 00) 


5960 CONTINUE 


H-*5GFUtJ>* 

(CDEL (EC(I,J)*ALB)+EC(I*U) ) ♦ ALB*TRM2/ 
((] *-AL0)* TAN(UtALB) ) ) 

SGC ( I • J) - SOMX'J) ♦ 8*»3GF(I.J>* 
(•TDEHECd. J) ,ALB) ♦EC(1*U1) * AL«*T«M2/ 
{ 1 1 • ♦ALBl * T ANH (U (ALB) ) 1 


C 


5975 CONTINUE 
5977 CONTINUE 
C 

C STRESS-STRAIN OUTPUT FOR THIS LOADING 

CALL OUT PS S (NORP tNTfl » SGO* SGC *£ V ) 

IFUNC*EO*l) GO TO t>060 
C RETAIN NUMBER OF FAILED MEMBERS 

IRT * ITOT 
C 

C CHECK COMPUTED STRESS AGAINST MAXIMUM ALLOWABLE STRESS 

6000 CONTINUE 
K ■ 0 

DO 6050 I » 1 • NOPP 
DO 6050 JMtNTR 
IF ( IFAIL ( I * J) > GO TO 6050 
IF ( ABS ( SGC < I « J) ) *LT • SGM < 1 1 J) ) GO TO 6050 
C STORE TEMPORARILY MEMBERS FAILING UNDER APPPlIED LOAD INCREMENT 

K • K 4 ] 

TMPFL ( I t JtK) » ABS (SGC < 1 • J) > 

U(K) « I 
JJ(K) • J 
6050 CONTINUE 

IF <K.EO*0> GO TO 6051 
C FIND FIRST FAILED MEMBER (OR MEMBERS) 

CALL RESET (TMPFLtKf I FA IL *SGM* 1 1 . JJtKMAA , I TOT,MEHB » 1 • OOE-02 • 

1 MXCAL * I SENT ) 

IF ( ISENT «EQ • 1 ) GO TO 7000 

C RETAIN SLOPE TERMS (IN TERMS OF FIRST FAILED MEM0EH) 

C FOR INTERPOLATION ROUTINE 


£> . 
A? A? 

•S 




$ 


,c? V 

4r 
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002120 

002122 

002124 

002132 


002136 

002136 

002146 

002155 

002155 

002156 

002160 

002162 

002164 

002164 

002166 

002167 

002172 

002176 

002203 

002205 


002210 


002217 

002217 

002221 

002223 

002225 

002231 

002233 


002233 

002233 

002235 

002240 

002241 

002242 

00224? 

002244 

002246 

002250 


!M»iI(KMAX> 

JM*JJ(KMaX> 

5GDEL«ABS (5GC ( IM* JH> -STORS ( IM* JM) 1 
SMDEL*$GM < IM* JM) -ABS (5T0RS ( JM* JM) ) 

C 

C CHECK SENTINELS GOVERNING FAILURE ANALYSIS 

6051 CONTINUE 

IF< <IT0T,EG.1RT> .AND. (MEMB.EQ.Z) > GO TO 6052 
C • T * — = NO FURTHER FAILURES AFTER FIRST REE VALUATION 

IFUNOT. < (ITOT.EQ.IRT) .AND* (MEMB.EQ.3) )) GO TO 6060 
C .F. — r NO FURTHER FAILURES AFTER MORE THAN ONE REEVALUATION 

6052 CONTINUE 
mEmb « 1 

DO 6055 1*1*6 

SLN ( I 9 1 1 * SGO(Itl) 

6055 CONTINUE 
6060 CONTINUE 

C STORE STRESS VALUES FOR PRESENT LOADING 
DO 6150 I»l*NORP 
DO 6150 J«1#NTR 
IF ( IFAll f I * J> ) 00 TO 6150 
STORS(IfJ) * SGC ( I * J) 

6150 CONTINUE 
C 

IF<INC*EQ.l) GO TO 6500 
C CHECK IF ALL MEMBERS HAVE FAILED 
IF ( ITOT ,GE.KNT) GO TO 7000 
C CHECK FOR FAILED MEMBERS 

C MEMB * 1 — «■ = NO FAILUHE OR ADDITIONAL FAILURE* INCREMENT IN 

C STANOARD FASHION 

C MEMB » 2 = FIRST FAILURE* CALCULATE S3* * REEVALUATE 

C MEMB ■ 3 —£ FURTHER FAILURES UNDER GIVEN SG*t CONTINUE REEVALUATION 

GO TO (6500*6200*015)* MEMB 
C 

C CALCULATE 5G*, INTERPOLATED LOAD 

6200 CONTINUE 

IF(INC.EQ.l) GO TO 6500 

DO 6250 1*1*6 

SLOEL * SLN (1*1) -SLP (1*1) 

SGO (1*1) » (SMDEL/SGOEL) # SLDEL ♦ SLP(Itl) 

6250 CONTINUE 
GO TO 915 
C 

C INCREMENT APPLIED LOAD 

6500 CONTINUE 

DO 6550 1*1*6 

SGO (1*1) * SGO t 1*1) * SM(I) 

6550 CONTINUE 
GO TO 915 
7000 CONTINUE 

C CONVERT RADIANS TO DEGREES, 

DO 219 I»l*NOR 
AL(I )*AL(I) *57,29578 
219 CONTINUE 

00 319 I*l*NORP 


vi 

Cn 



Ln 

<Xi 
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002252 DO 319 Jwl*NTR 

002253 BT 1 1 ♦ J> «BT < 1 t J> *57 *29578 

002257 319 CONTINUE 

002263 60 TO 1 

C 

c 

c 

002264 9710 FORMAT <//* MATRIX IS SINGULAR OF RANK ■ **12, 

1 *+ DETERMINANT * **E15*5/> 

002264 9715 FORMAT <6E12«5) 

002264 9720 FORMAT <//* SOLUTION FOR STRESS CONVERGES In «tI2* 

1 * ITERATIONS*//) 

C 

c 

c 

c 

002264 33 STOP 

002266 END 
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ooooio 

ooooio 

ooooio 


OOOOIO 

000015 

000021 

000025 

000031 

0000^3 

000051 

000064 

000072 

000076 

000103 

000104 

000105 

000157 

000224 

000231 

000235 

000241 

000245 

000252 

000253 

000334 

000410 

000415 

000421 

000425 

000425 

000425 

000425 

000425 

000425 

000425 

000425 


000425 

000425 


C 


C 

C 


SUBROUTINE INPUT <NOR *NTP » ICS* I Af_« IBT) 

COMMON /C INPUT/ AD»BD»DD*NB*L *EY tEC i A * VHOD* VTI • V* SON 

REAL NB( 5i 30 > *L ( 3»30)tl8T( 5*3o>*IAU 5) 

DIMENSION EY ( 5*30>*A< 5f30)«EC( 5»30U 

1 VTH 5*30)*SGM< 5t 30 > 


WRITE (6*9610) ICS 
WRITE (6*9615) 

WRlTE<6*9620) 

WRITE (6*9615) 

WRITE (6*9625) A0*BD*DD 
WRITE (6*9627) V 
WRITE £6*9630) NQR*NTR 
WRITE (6*9631) VMOO 
WRITE (6*9635 ) 

NORP * NOR ♦ 1 

DO 100 I -ltNORP 

DO 100 J «1 *NTR 

IF 1 1 .NE.NORP) WRITE ( 6*9650) 

1 

IFd.EO.NORP> WRITE (6*9651 1 

1 

100 CONTINUE 

WRITE (6*9640) 

WRITE (6*9615) 

WRITE <6*9645) 

00 150 I wltNORP 

DO 150 J -ltNTR 

IFd, NE.NORP) WRITE 46*9650) 

1 

lFd.EO.NORP) WRITE [6*9651 ) 

1 

150 CONTINUE 

WRITE (6*9660 ) 

WRITE <6.9615) 


IALd> tieidf J)fNB(ItJ) «L(Iivl)t 
VTHItJ) 

IBT d.J) *NB ( I » J) «L ( I * J) t 
VTld.J) 


lALd) tIBTd»J)tEY(ItJ) ,EC(ItJ) * 
A(IiJ) *SOM(I, J) 

IBT d«J)fEYd«J) «ECdt J» • 
AII.J) *SG«d*J) 


c 

9610 FORMAT 
9615 FORMAT 
9620 FORMAT 

9625 format 

1 

962T FORMAT 

9630 FORMAT 

1 

9631 FORMAT 
9645 FORMAT 

1 

2 
3 

9640 FORMAT 
9635 FORMAT 
1 


<1H1 i50X* # SPACE FRAME MODEL** 12 > 

( 51X ** -*) 

<//5lX f *MdOEL GEOMtTRY*) 

</20***DlMENS10NS OF MODEL** 1 0* **A "**E12. 5*5X**ft «**E12.5t 
5X»*C ■**£12.5) 

<5oX**VOLUME ■* tE | 2.5) 

</2OX**NUM0ER Of I KUSS ORIENTATIONS **#I2i 
/2OX.*NUM0ER OF FIBER ORIENTATIONS P£R TRUSS »**I2) 
(/20X # *VOLUME FRACTION Of MODEL «**E12.5) 

(/29X»*0RIENTAT10N5* *20* *• YOUNG! S* * 1 OX **ECCENT •* * 

10X.«CROSS-**10Xi*MAXIMUM*/26^,*TRUSS*tlOX**FlBERS IN** 
lOXt*MOOULUS**25X,*SECTIONAL*f 12X**STRESS*/43X**TMUSS** 
47X**APEA*) 

< //50X * *F IBER PROPERTIES*) 

</29X f *0RIENTATI0NS**18X**NUM8tH OF**TX* *LENGTH OF*t 
aX**VOLUME** 
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INPUT 
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2 /26Xf*TRDSS**10X**FlREH5 IN* t 13X'*6AYS*t 13* » •TRUSS** 

3 1 0* »*FRACT ION** 

4 /45X,*TRUSS*/) 

000425 9650 FORMAT ( 24X *2 <£12,S*2X 3 ,4 < 5X*E12.5» 

000425 9651 FORMAT ( 38X tEl 2.5*2X t 4( 5XfE12*5» 

000425 9680 FORMAT (//50X**ANAlYSIS*) 

000425 RETURN 

000426 END 
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000017 

000017 

000017 

000017 

000017 

000020 

000022 

000024 

000026 

000037 

000041 


000044 

0000$1 

000057 

OOOOM 

000063 

000066 

000071 

000072 

000112 


000120 

000123 

000125 

000130 

000132 

000137 

000142 

000144 

000144 

000144 


SUBROUTINE RESET <TMPFL*K,IFAIL*SGMt II *JJ*KMAX*IT*MtMe f £PS.MXCAL* 

* ISENT ) 

C0MM0N/L0C2/ STO 
C 

LOGICAL IFAILC 5*301 

PEAL RAT (100) *TMPFL( 5*30*40 ) *SGM ( 5* 30 ) 9 L ARGE*EPS 

INTEGER II (40) • JJ(40) *510(5*30) 

ISENT * 0 

C FORM RATIO OF STRESS AT FAILURE TO MAXIMUM STRESS 
00 50 1*1 *K 
iM.tim 

RAT ( I ) »TMPFL ( IM* JM*I ) /SGM ( IM, JM) 

50 CONTINUE 

C FIND LARGEST RATIO 

CALL MAXVAL(K'RATtKMAX'LARGF) 

C DETERMINE FAILED MEMBERS FOR THIS LOAD 

C (I.E* THOSE WITHIN EPS OF MAXIMUM RATIO) 

DO 100 I-ltK 

IF ( (LARGE-RAT (I)) *GT,EPS) GO TO 100 

itemp ■ inn 

JTEMP • JJ(I> 

IFAIL(ITEMPtJTEMP) * .TRUE* 

STO ( ITEMP. JTEHP)-! 

IT « IT ♦ 1 

WRITE (6*9720) II(I)*JJ(I> 

100 CONTINUE 
C 

C SET REEVALUATION COUNTER ■ 0 BEFORE INITIAL REEVALUATION 
1FIMEMB.EQ.2) KRCAL°0 

C RESET SENTINEL ONLY AFTER FIRST REEVALUATION 
IF (MEMB.EG* 2> MEMB*3 

C SET SENTINEL INDICATING FIRST MEMBER FAILURE UNDER PRESENT LOAD 
IF (MEMB*EO» 1 ) MEMB«2 
IF (MEMB.NE*3) RETURN 

C LIMIT NUMBER OF REEVALUATIONS PE« LOADING BY *MXCAL* 

IF (KRCAL^GE.MXCAL) isent-i 
IF (KRCAL.GE.MXCAL) RETURN 
KRCAL ■ KRCAL ♦ 1 
RETURN 

9720 FORMAT ( 1 OX**M£MBER FAILURE 1 **3X** TRUSS* *I2*3X**F I HER* *12) 

END 
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000010 

000010 

0000)0 

0000)0 

oooolo 

0000)0 

000013 

U00020 

□00021 

00002* 

000026 

00003? 

000034 

000035 

000043 

000045 

000071 

000076 

000106 

000117 

000121 

0001 ?? 

000124 

000125 

000147 

000205 

000212 

000212 

000212 

000212 


000212 

000212 

000212 


SUBROUTINE OUTPSS (NORP * NTR* SGO »SSF *£V J 
C 

COMMON/L.OC2/ STO 

REAL SGO <6tl)fEV(6'l)t SOP l 5*30) 

DIMENSION DR (6) 

INTEGER STO (5»30 ) 

C 


DAT A DR/4H 11, 4H 22t4H 33, 4H 12*4H 23>4H 31/ 

WRITE (6,9750 ) 

C AFTER A FIBER FAILS SET THE STRESS TO ZERO, 

00 2 K J3*l ,5 
00 2 KJ4*1,30 

IF ( STO ( K J3 , K J4 ) • EG, 0) GO TO 2 
SGF(KJ3*KJ4>*0*0 
2 CONTINUE 

C SET ALL ELEMENTS OF MATRIX STO TO ZERO. 

DO 1 KJ1-1*5 
00 1 KJ2«1,30 
1 STO<KJ1,KJ2>*0 
DO 100 1*1,6 

WRITEi 6, 9762) DR < I ) , SGO I J , 1 ) ,EV< 1,1) 

100 CONTINUE 

IF (NTR »EO. 1 > WR1 TE < 6 » 9769) 

IF(NTR.GT.l) WRITE <6*9770 > 

DO 150 1*1 »NQRP 
DO ISO J*1*NTR,2 
JR * J*1 

NMJ » NTR-J 

IF (NHJ*EG« 0 ) WRI T6 < 6 #9775) I,Ji5GF<I,J> 

IF<NMJ.GT*0> WRITE (6,9775) I , JfSGF <| , J) , I , JP,SGF < I « JP) 
150 CONTINUE 


9750 FORMAT 
9762 FORMAT 

9769 FORMAT 
1 

9770 FORMAT 
1 

2 

9775 FORMAT 
RETURN 
END 


< / /20X **DIRECT ION*t 20X ** APPLIED loAD*,20Xi *STRA IN*/ ) 
(?4X,A4,2<20XtEl5,5> > 

(?OX,*ORlENTATIONS^tl5X,*STRESS*/10X»*T«USS*,iOX, 
♦FIBERS IN*/27X,*TRUSS»/1 

l?(2OX,*ORlENTATlONS*,I5Xf*STRESS*t5X)/l0X.«TRUS5**l0X» 
•FIBERS IN* 1 30X * • TRUSS*, 1 QX t*F IRERS IN*/2/X f #TRUSS*i 
50X,*TRUSS*/) 

(?U2X,I2ti2X,lZ*l5X,E15.5) ) 
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OOOOlt 

U00Q11 

000011 

UOOOll 

000013 

000030 

000034 

0000*0 

000053 

000060 

OOOO^l 

UOOllO 

000120 

000120 

000120 

000120 

000120 

000120 

000120 


SUBROUTINE COUMP(NORP*NTR*ICON*IT*SGF*SGP> 

C 

REAL S6F< 5*30)*SGP( 5*303 

0IMEN5I0N C0NMSe<2*2» 

DATA CONMSG/lOH DOES NOT *10HCONVERGE *10H DIVERGES » 

« 10H / 

ITM » IT-1 

WRITE <6 *9710) (CONMSG ( Jt ICON! * Jm\ .4) 

WRITE ( 6*9720 ) 

WRITE (6*9725) 

WRITE (6*9726) IT*ITM 
DO 100 Ul*NORP 

DO 100 J«1*NTR 

WRITE (6*9730) I ♦ J* SGF ( I * J) *SGPCI*J> 

100 CONTINUE 

9710 FORMAT t//l 0X*#$OLUTION FOR STRESS *,4A4//> 

9720 FORMAT ( 1 OX**TERNINAL STRESS VALUES*) 

9725 FORMAT { 10X**TRUSS** 1 0* * • FIBERS *i 10* * 2( 1 OX,#STRESS*> ) 

9726 FORMAT (56X*2 (* ( ITER*** 12**) ** 10X) ) 

9730 FORMAT ( 12X* 12* 1 5X * 1 2*l5X *2 (5X*E1 2*5) ) 

RETURN 

END 



cr\ 

to 
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000007 

000007 

000007 

000010 

000012 

000013 

000016 

000017 

000020 

000023 

000023 

000024 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


SUBROUTINE MAXVAL<K#Y* KMAX, UPPER) 

DETERMINE maximum of array OF VALUES 

AND SUBSCRIPT ASSOCIATED XI TM MAXIMUM 

PARAMETERS USED 


Y I ARRAY NAME 

K t NUMBER OF ELEMENTS IN ARRAY 

UPPER I MAXIMUM VALUE OF Y 

KMAX 1 SUBSCRIPT ASSOCIATED WITH UPPER 


REAL Y ( K ) » UP PER 

UPPER » Y(l) 

KMAX ■ 1 

IF (KtEQ»l ) GO TO 150 
DO 100 1 b2 »K 

IF (Y (1 > ♦LT* UPPER) 00 TO 100 
UPPER « Ytn 

KMAX « I 
100 CONTINUE 
150 CONTINUE 
RETURN 
END 
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OOOOU 

C 

SUBROUTINE rtXttULO(A'B'C'NNOWAtNCOLA'NCOLB'MA,NA*NB> 
DIMENSION A(NtiOWAfNCQLA) (NCOLA tNCOLB) ,C (NRO*A tNCULB) 

RQ6DOOO? 

000014 


REAL A*B.C*X 

R0600003 

000014 


DO 20 I-l.MA 

R06D0004 

000015 


DO 20 uMltNB 

RU600005 

000016 


x«o. 

RC6D0006 

000017 


DO 10 K»ltNA 

R06D0007 

000021 

10 

X»X*A(I,K)»BtK* J) 

R06D0CCB 

000035 

20 

C(1*J)-X 

R06D0009 

000047 


RETURN 


000047 


END 
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000012 

000032 

000012 

000013 

00001* 

000027 

00003* 

000035 

000036 

000037 

000040 

000052 

00005* 

000041 

000063 

00006 * 

000066 

U0006T 

000072 

000073 

00007* 

000075 

000076 

000077 

oonioo 

000102 

000103 

000105 

000112 

ooon* 

000116 

000H7 

000120 

000123 

000126 


SUBROUTINE INVRTD{A*NDIM*N*PETA*EP$* 1 RANK » UNDER ) 

c 

C ROUTINE INVHTD INVERTS AN N X N MATRIX USING GAUSS-JORQAN 

C ELIMINATION METHOD 

C 

c * VARIABLE DICTIONARY * 

c 

C A(NiN) 

C NOIM 

C N 

C TRANK 

C DET A 

C EPS 

C 

C UNDER 

C 

c «M*i# 

c 


I MATRIX TO BE INVERTED PASSED * INVERSE RETURNED 
1 UPPER LIMIT TO MATRIX DIMENSION 
I DIMENSION OF MATRIX 
I RANK OF MATRIX 

I deterinant of matrix 

I ADJUSTABLE TOLERANCE FACTOR COMPARED TO 
VALUE OF PIVOTAL element uuring inversion 
1 UNDERFLOW LIMIT (CHECK ON COMPUTED VAR*) 



DIMENSION A(NDIM*NDIM> 





INTEGER IR(60)«IC(60)»R»S 





CHECK MATRIX ELEMENTS FOR 

UNDERFLOW 

POSSIBILITIES 



DO 5 1*1 





DO 5 JaJ »N 





IF C ABS ( A ( I * J) ) *LT* UNDER) 

A(IfJ) ■ 

0*0E 00 


5 

CONTINUE 





DETA-1. 



R12D0006 


SUMaO* 



R12D0007 


DO 10 I * 1 «N 



R12D0008 


no io J « l § n 



R1200009 

10 

SUM«SUM* A ( I » J) **2 



R12D0010 


SUM. SORT (SUM) 



R12D001 1 


DMA • N**Z 



R12D0012 


PMS*5UM/DMA 



R12D0013 


tol»eps*rms 



R12D0014 


DO 20 I • ItN 



R1200015 


IR(I)*0 



R1200016 

20 

ICm«0 



RI2D0017 


S«0 



R1200018 


R * N 



R 1200019 

30 

I“0 



R12D0020 


J*C 



R12000Z1 


TESTaO.O 



R12D0022 


DO 50 K * 1»N 



R1200023 


IF ( IRIK) ,NE*0)GO T050 



R1 20002* 


no *o L * i • n 



R 120 0025 


IF(IC(L) *NE*0)GO TOAO 



R1200026 


X* ABS (A (Kf L) ) 



R12D0027 


IF(X*LT.TEST)G0 to*o 



R12O0028 


1 *K 



R1200029 


J«L 



R1 200030 


test»x 



R12D0031 

*0 

continue 



R1200032 

50 

CONTINUE 



RI2D0033 


PI V*A ( I * J) 



R12D003* 


sg 

00 

88 


.O 

%/ 





RUN VERSION g.3 — PSR LEVEL 363 


INVRTO 


000133 IF C ABSfDETA) .LT. UNDER) DETA * O.OOE 00 

000137 deta«piv#deta 

0001*1 IF ( ABSCPIV) «LC. TOL) GO TO 150 

0001*3 IR (I > »J 

0001*5 IC<J)«I 

0001*6 PlV ■ l.OEO/PIV 

000150 All*J>»PtV 

000153 00 60 K m liN 

000155 60 IFIK.NE.J) A<I»K)«A(I*K)*PIV 

000166 DO 90 K * liN 

000167 IF (K.EQ.I) GO TO 90 

000171 PIV1 • A (K * J) 

000175 70 DO 80 L * l'N 

000177 IF { A0S<P1V1) # LT. UNDER) PIV1«0,0E 00 

000203 IF { ABS<A<I*L) > .LT, UNDER) A<I*D « 0*00E OD 

000217 SO IFCL.NE. J) A<K»L) B A(K*U-PIVl*A(If L> 

000235 90 CONTINUE 

0002*0 DO 100 K • 1 * N 

0002*1 100 IFCK.NE.I) A<K, J)*-PlV*A<Kf J) 

000253 S*S*1 

00025* IF(S.LT.R)GO TO 30 

000256 110 DO 1*0 X » liN 

000260 K«ICm 

000262 R*lRtn 

000263 IF ( K.EQ. I ) GO TO 1*0 

000265 DETA--DETA 

000266 DO 120 L « IfN 

000267 TEMP-A £K,L> 

000273 A(K.L)*A(1»L) 

000302 120 A ( I tL ) *TEMP 

000306 DO 130 L ■ ltN 

000307 TEHP-A (LtM) 

00031* A<Lt*>»A(Lf I) 

000322 130 A (L 1 1 ) *TFMP 

000330 IC(M)»K 

000332 IR<K>»M 

000333 1*0 CONTINUE 

000336 150 IRANK«S 

000337 RETURN 

0003*0 END 


ip 
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Rl 2D0035 
R1 200036 
R12D0037 
R1200038 
R 12D0039 
R 12000*0 
R12DOO*! 
P12D00*2 
R12D00*3 
R12D00** 
R12DOO*5 
R12DOO*6 


R12DOOA7 
Rl2D00*6 
Ri2DOO*9 
P12D0050 
R12D0051 
R1200052 
R1200053 
R12D005* 
R12D0055 
R12O0056 
R12O0057 
R12DOO50 
R12D0059 
R12D0060 
R1200061 
Rl 200062 
R12D0063 
R12DQ06* 
R12D0065 
R12D0066 
R12D0067 
R12D006B 
R12D0069 
R12D0070 
R12D0071 





